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THE  DEVELOPMENT  OF  TECHNIQUES  FOR  OPERATIONAL  EXPLOITATION  OF 
SATELLITE-OBSERVED  SOUNDING  DATA  FOR  RESOLUTION  OF  ATMOSPHERIC 
THERMAL-STRUCTURE  VARIABILITIES 

Phase  Two:  Techniques  for  Processing  VTPR  Side-Scan  Radiances  Including 
Diagnosis  for  the  Clear-Column  Radiance  Components 


1  •  Summary  of  Developments 

1 . 1  Developments  of  Phase  One 

The  major  developments  performed  under  the  earlier  Phase  One 
include: 

(1)  A  generalization  of  the  mass-structure  model  to  any  specified  selection 
of  standard  pressure  levels.  The  mass-structure  model  constitutes  the 
information  framework  for  upper-air  analysis,  including  resolution  of 
the  atmospheric  thermal  structure. 

The  degrees  of  freedom  of  the  mass -structure  model  are  defined  as 

follows:  The  atmospheric  column  is  divided  into  N  layers,  or  segments, 

with  pressure  as  vertical  coordinate.  In  each  segment  the  virtual 

temperature  is  linear  in  p*,  where  k  =  R/c  .  This  makes  the  chosen 

P 

form  of  static  stability  constant  in  each  layer,  changing  abruptly  at  the 
interfaces.  The  structure  model  incorporates  the  hydrostatic  and  gas-law 
approximations.  With  pressure  as  vertical  coordinate,  the  relationship 
betv/een  heights ,  temperatures,  and  layer  stabilities,  is  linear.  This 
linearity  is  the  virtue  of  the  mass -structure  model.  The  transformation 
matrices  are  produced  by  our  generalized  program. 

The  generalized  program  permits  arbitrary  selection  of  a  set  of  N+2 
pressure  levels,  referred  to  as  the  standard  pressure  levels,  in 
monotonic  order.  The  program  locates  the  interfaces  midway  in  pressure 
between  the  standard  levels  except  for  omitting  interfaces  between  the 
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first  two  and  the  last  two.  This  defines  N  layers,  for  the  desired 
vertical-resolution  powers. 

A  ten-level,  eight-layer  model,  based  on  the  traditional  set  of 
standard  pressure  levels,  has  been  in  operational  use  at  the  Fleet 
Numerical  Weather  Central  (rNWC)  for  many  years.  The  resolution  powers 
require  supplementation  in  the  1000-800  and  100-0  mb  ranges. 

(2)  Radiances  transformed  into  terms  of  the  modelled  thermal  structure.  The 
production  of  transformations  which  express  each  channel  radiance 
anomaly  in  terms  of  a  linear  combination  of  the  atmospheric  thermal- 
structure  anomalies  including  a  sea-surface,  or  ground,  temperature 
anomaly. 

An  arbitrary  standard  thermal  structure  is  defined,  including  a  sea 
surface  or  ground  temperature.  Given  a  transmission  function,  the 
standard  radiance  is  calculated  for  that  channel  in  terms  of  an 
atmospheric  contribution  and  a  ground  contribution.  This  is  done  by 
integrating  the  radiative  transfer  equation  for  the  defined  standard 
atmosphere. 

The  thermal  structure,  being  linear-in-p*  in  each  layer,  now  has 
its  degrees  of  freedom  defined  in  terms  of  the  temperature  anomaly  at 
each  interface  and  at  the  base  (1000  mb)  and  top.  A  ground 
temperature  anomaly  is  also  permitted.  The  Planck  function  is  linearized 
about  the  assumed  standard  temperature  structure,  enabling  transformation 
of  the  radiative-transfer  equation:  the  radiance  anomaly  is  defined  as 
a  linear  combination  of  the  temperature  anomalies. 

The  mass-structure  transforms,  produced  by  program  (1)  above, 
enable  the  temperature  anomalies  to  be  expressed  in  terms  of  other  sets 
defining  the  therma  1-structure  anomaly.  This  permits  expressing  the 
radiance  anomaly  in  alternate  forms:  as  a  linear  combination  of  height 
anomalies  of  the  standard  pressure  levels,  as  a  linear  combination  of 
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layer  thicknesses,  as  a  linear  combination  of  layer  stabilities  and  a 
thickness  (e.g.,  1000-300  mb),  and  in  other  forms. 

For  a  specified  transmission  function  and  an  assumed  standard 
atmosphere,  this  program  generates  the  standard  radiance  contributions 
from  atmosphere  and  ground,  and  several  sets  of  combination  coefficients 
which  express  the  radiance  anomaly  in  terms  of  corresponding  sets  of 
thermal-structure  parameter  anomalies. 

These  transforms  apply  to  clear-column  radiances,  as  do  the 
specified  transmission  functions  on  which  they  are  based. 

(3)  A  generalized  retrieval  scheme  for  transforming  clear-column  radiance 
estimates  into  enhanced  resolution  of  the  atmospheric  thermal  structure 
and  the  ground,  or  sea  surface,  temperature. 

This  retrieval  is  in  the  form  of  a  structure  blending  scheme.  The 
input  information,  expressed  by  weighted  clear-column-radiance 
estimates  and  by  independent  weighted  estimates  of  a  set  of  atmospheric 
thermal-structure  parameters,  is  blended  to  produce  the  optimum 
resultant  implications,  with  associated  resolutions  afforded  by  the 
input  information.  The  weights  are  defined  as  the  inverse  of  the 
associated  error  variances. 

Incorporated  into  this  blending  scheme  is  the  absolute  constraint 
that  all  static-stability  implications  be  positive.  This  constraint 
represents  additional  input  information. 

The  information  blending  methodology,  in  general,  combines 
independent  weighted  estimates  into  the  non-independent  resultants 
implied  by  the  focus  of  all  the  input  information.  It  handles  linear 
interrelationships  and  other  extensions.  The  resolution  that  comes 
out  reflects  what  wont  in.  In  the  present  application  of  the  methodology 
we  output  not  only  the  resultant  implication  for  each  thermal-structure 
parameter  but  also  for  each  radiance. 
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The  blending  can  be  forced  to  perfectly  accommodate  fill  input 
radiance  estimates*  by  setting  the  associated  input  weights  very  high 
compared  to  any  conflicting  information  in  the  independent  thermal- 
structure  estimates.  But  realistic  weighting  is  much  to  be  preferred. 

The  blending  scheme  includes  three  cycles  of  retrieval  in  order 
to  enable  roevaluation  of  the  input  estimates  and  input  weights  in 
comparison  with  resultant  indirect  implications.  Estimates  judged  to 
be  in  gross  error  are  rejected,  and  the  weights  of  maverick  estimates 
are  reduced.  The  constraint  that  the  retrieved  static  stabilities  be 
positive  amplifies  the  disparity  with  contaminated  estimates.  The 
nature  of  the  structure-blending  scheme  exploits  the  radiance 
information  for  adjustment  of  the  related  structure  components  in 
conformance  with  the  firmness,  and  lack  of  firmness,  in  the  input 
estimates.  That  is,  the  transformation  of  the  input  information  is 
dictated  by  the  specific  set  of  input  weights,  and  by  the  analysis 
framework  and  the  interrelationships.  Various  examples  of  this 
interplay  have  been  demonstrated. 

The  Phase  One  research  effort  is  documented  in  the  following  Task 
Performance  reports: 

Task  One:  "The  development  of  techniques  for  operational  exploitation 

of  satellite  infrared  channel  radiances  in  the  resolution  of 
atmospheric  thermal-structure  variabilities — task  one: 
perspective",  Technical  Report  of  Task  One,  Project  to. -181, 
Contract  No.  N62306-72-C-0125,  September  1972, 
Meteorology  International  Incorporated ,  Monterey,  California, 

61  pp. 

Contents:  Nature  of  the  Radiances  and  the  Retrieval  Problem, 
Context  for  the  use  of  Radiance  Information,  Review 
of  Retrieval  Schemes  with  Comments. 


Provided  that  these  estimates  do  not  inherently  violate  the  static-stability 
constraint. 


Task  Two:  "The  development  of  techniques  for  operational  exploitation 

of  satellite  infrared  channel  radiances  in  the  resolution  of 
atmospheric  thermal-structure  variabilities--ta.sk  two:  new 
formulations".  Technical  Report  of  Task  Two,  Project  M-181, 
Contract  No.  N62306-72  -C-01 25  ,  March  1973,  Meteorology 
International  Incorporated ,  Monterey,  California,  91  pp. 

Contents:  Mass-Structure  Model,  Parameterization  of  Thermal 
Structure  in  Terms  of  Layer  Static  Stabilities, 

Transformation  of  Radiances  in  Terms  of  Linear 
Combinations  of  Thermal-Structure  Parameters,  Retrieval 
by  Structure  Blending. 

Task  Three:  "The  development  of  techniques  for  operational  exploitation 

of  satellite  observed  sounding  data  for  resolution  of  atmospheric 
thermal-structure  variabilities — task  three  and  final  report". 
Technical  Report  of  Task  Three,  Project  M-181,  Contract  No. 

N 62306-72-C-01 25  ,  June  1973,  Meteorology  International 
Incorporated,  Monterey,  California,  59  pp. 

Contents:  Generalization  of  the  Mass-Structure  Model,  Corresponding 
Transformation  of  Radiances  in  Terms  of  Linear  Combinations 
of  Thermal-Structure  Parameters,  Generalization  of  Structure 
Blending  Including  Constrained  Blending  (Assuring  a 
Statically-Stable  Retrieval)  and  Reevaluations  of  Estimates 
and  Weights . 
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1 . 2  Developments  of  Phase  Two 


The  major  developments  performed  under  Phase  Two,  and  detailed  in 
the  present  Technical  Summary  Report,  include: 

(1)  Generalization  of  the  mass-structure  model  and  the  transforms  for 
relating  clear-column  radiances  to  the  atmospheric  thermal  structure. 
This  program  represents  a  merging  of  items  (1)  and  (2)  listed  under 
Phase  One,  above,  but  with  added  generalizations  and  extensions. 

A  general  capability  has  been  developed  which  produces  standard- 
radiance  values,  and  the  transforms  which  express  radiance  anomalies 
in  terms  of  linear  combinations  of  atmospheric  thermal-structure- 
parameter  anomalies.  This  capability  is  general  in  that  it  permits 
arbitrary  specification  of  (a)  a  table  of  standard  pressure  levels  which 
define  the  degrees  of  freedom  of  the  thermal-structure  model,  (b)  a 
standard  temperature  structure  in  terms  of  values  for  the  degrees  of 
freedom  of  the  modelled  thermal  structure,  and  (c)  the  transmission 
function,  for  any  channel  and  nadir  angle,  tabled  in  terms  of  a  specified 
increment  in  p*. 

This  general  capability  was  applied  to  the  transmission  functions 
of  an  operational  Air  Force  VTPR  sounder.  We  were  provided  with  the 
transmission  functions  for  a  set  of  channels  and  a  range  of  nadir  angles. 
The  general  program,  and  the  tables  produced  for  an  18-level  structure 
model  and  a  specified  standard  thermal  structure,  are  included  in  the 
deliverable  items.  The  tables  of  transforms  were  used  in  several 
contexts  including  the  normalization  of  radiances  for  any  nadir  angle. 

We  also  investigated  the  designed  differential  resolution  powers  of  the 
sounder  afforded  by  the  specifications  of  the  channels  over  the  range 


(2)  Examinations  of  real  satellite  radiances,  and  investigation  of  the 
characteristics  of  cloud  contamination. 

We  were  also  provided  with  a  sample  of  radiances  measured  by  the 
Air  Force  VTPR  sounder,  taken  in  orbit  on  26  March  1974.  All  complete 
side-to-side  scans  on  this  sample  tape  were  transformed  into  normalized 
radiance  anomalies  for  each  channel  and  nadir  angle.  Each  scon  was 
plotted  (by  Varian  plotter)  in  the  form  of  a  cross  section,  placing  each 
normalized  anomaly  at  the  center  of  its  energy  source  region.  These 
cross  sections  were  studied  in  our  investigation  of  the  character  of 
cloud  contamination. 

(3)  Development  of  a  capability  for  diagnosing  the  clear-column  radiance 
components  (Program  CLRX) . 

The  study  of  radiance  cross  sections  inspired  the  conceptions  and 
the  major  achievement  of  Phase  Two:  a  capability  to  diagnose  each 
complete  scan  for  the  inherent  resolution  of  the  clear-column  radiance 
component  of  each  measured  radiance.  We  refer  to  this  programmed 
capability  as  CLRX.  The  program  operates  on  the  scan  after  it  has 
been  processed  by  the  preceding  development,  item  (2)  above;  however, 
it  can  be  modified  to  operate  directly  on  the  measured  radiances. 

CLRX  produces  not  only  an  estimate  of  the  clear-column  component 
of  each  radiance  element  of  the  scan  but  also  an  associated  reliability 
for  each  such  estimate. 

During  the  latter  part  of  the  performance  period  we  worked  on 
producing  a  comprehensive  version  of  CLRX  for  use  in  oceanic  regions. 

In  order  to  demonstrate  and  assess  the  potential  powers  of  CLRX  we  did 
not  input  any  independent  thermal-structure  information  in  the  form  of 
first-guess,  or  background,  estimates  of  the  clear-column  radiances. 

We  wished  to  demonstrate  the  powers  of  CLRX  to  diagnose  the  information 
inherent  only  in  the  measured  radiances  of  the  scan. 
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CLRX  represents  a  significant  breakthrough  in  the  practical 
exploitation  of  satellite  V'i’PR  sounders.  CLRX  is  not  dependent  on 
independent  sea-surface~temperature  (SST)  values  for  categorical 
calibration  of  the  clcar-column  radiances.  Hence,  its  present 
applicability  extends  into  all  oceanic  regions,  even  where  there  is 
little  independent  SST  information.  Wo  have,  in  fact,  conceived  of 
the  extension  of  CLRX  for  use  over  land  regions.  The  formulations 
promise  to  re  much  more  complex  even  for  terrain  of  low  and  uniform 
elevation.  Computer  running  time,  however,  would  be  comparable. 
Extension  for  application  in  regions  of  high  and  irregular  elevations 
is  also  conceivable  for  diagnosis  of  the  lesser  penetrating  channels. 

Processing  of  the  measured  radiances  by  wLRX,  for  diagnosis  of 
the  clear-column  components,  is  essentially  independent  of  the  accuracy 
of  the  specified  transmission  functions  for  the  channels  at  nadir  angles. 
We  have  not  unduly  involved  ourselves  with  the  accuracy  of  the 
transmission  functions,  with  adjustments  for  attenuation  by  water  vapor 
and  aerosols,  and  with  the  quality  control  and  calibration  of  channels. 
These  matters  have  been  outside  the  scope  of  our  project.  The 
techniques  that  we  have  developed  are  adaptable  to  specification  of 
transmission  functions  relative  to  individual  side-to-side  scans. 

Complex  refinements  of  the  transmission  functions  can  be  entertained. 

CLRX  produces  estimates  of  the  clear-column  radiance  components 
across  the  entire  scan,  thus  retaining  whatever  gradient  information  is 
inherent  in  the  measured  radiances  of  that  scan. 

CLRX  is  based  on  a  model  for  interpreting  cloud  contamination . 

The  largest  differences  that  arise  in  the  radiances  between  two  spots 
in  proximity,  are  generally  due  to  differences  in  the  amount  of  cloud 
at  only  one  cloud-top  level.  Based  on  this  interpretation,  large 
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differences  in  radiances  between  nearby  spots  of  a  scan*  produce 
estimates  of  the  local  correction  axis  (i.e.  ,  unit  vector)  for  adjusting 
the  measured  radiances  toward  their  dear-column  components.  CLRX 
does  not  require  the  cloud  top  to  be  at  the  same  level  over  the  entire 
scan;  a  gradual  slope,  or  even  stepwise  changes,  are  tolerated  in  the 
exploitation. 

It  should  be  understood  that  CLRX  can  only  exploit  information  that 
is  inherent  in  the  measured  radiances.  The  yield  of  information  is 
measured  by  the  reliability  estimate  which  is  produced  with  each 
diagnosed  clear-column  radiance  component.  The  yield  is  largely  a 
function  of  the  character  of  the  cloudiness  over  the  scan.  The  diagnosis 
is  generally  more  reliable  the  greater  are  the  contrasts  of  the  clear 
proportions  between  spots.  A  solid  cloud  top  over  the  entire  scan 
defies  all  diagnosis  of  the  affected  channels.  This  variability  in  yield 
underlines  how  essential  it  is  to  produce  the  associated  reliability 
measures . 

(4)  A  retrieval  demonstration. 

We  have  now  dealt  with  all  technical  impediments  to  the  operational 
exploitation  of  the  inherent  information  in  radiances  measured  by 
scanning  VTPR  satellite  sounders.  The  techniques  we  have  developed 
address  not  only  absolute  temperature  information  but  also  information 
for  resolution  of  relative  properties  such  as  horizontal  gradients.  Our 
techniques  are  not  transforms  based  on  idealized  theories.  They  are 


The  VTPR  scans  of  NOAA  spacecraft  are  more  closely  spaced  along  the 
orbital  path,  affording  useful  information  redundancies.  This  redundancy 
warrants  exploitation  by  analyzing  differences  between  measured  spots  not 
just  within  one  scan  but  also  between  adjoining  scans. 
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designed  for  use  in  the  real  world  of  the  operational  context,  and 
include  the  ability  to  cope  with  the  concomitant  problems  of 
contributing  normal  error  variances  and  gross,  anormal  errors  in  the 
data  streams. 

In  order  to  demonstrate  the  essential  completeness  of  the  range  of 
techniques  developed  under  Phases  One  and  Two  we  have  submitted 
output  from  CI..RX  into  the  generalized  structure-blending  retrieval 
scheme  developed  under  Phase  One.  This  retrieval  scheme  accepts 
weighted  clear-column  radiances,  and  transforms  this  information  into 
enhancement  of  the  resolution  of  the  thermal  structure  of  the  atmospheric 
column.  This  comprehensive  capability  also  exploits  the  information 
that  the  real  structure  is  statically  stable. 

This  demonstration  required  modification  of  the  structure-blending 
retrieval  scheme  from  the  13-lcvel,  1 1-layer  atmospheric-structure 
model  used  in  Phase  One  to  the  18-level,  16-layer  model  used  in 
Phase  Two. 

(5)  The  design  of  a  generalized  capability. 

A  test  and  evaluation  phase  is  now  essential  for  evaluating  the 
inherent  information,  in  VTPR  scanning  sounders,  that  can  be  exploited 
by  the  developed  techniques  under  operational  conditions.  We  have 
designed  the  generalized  basic  capability. 

This  capability  involves  interfacing  with  an  analysis  and  prediction 
system,  including  upper-air  thermal-structure  parameters,  and  analyses 
of  soa-surface-temperature  fields.  Provisions  include  the  collection 
of  high  quality  (i.e.  ,  of  the  highest  associated  weights)  CLRX-diagnosed 
clear-column  radiance  estimates;  values,  weights,  position  and  time 
data  would  be  saved  for  calibration  studies  with  any  available 
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independent  thermal-structure  information.  Supplemental  component 
could  be  systematically  incorporated.  Competitive  components  and 
variations  could  be  evaluated  side  by  side.  Realization  of  this 
capability  will  enhance  prospects  for  systematic  improvements  and 
extensions,  including  use  over  land  as  well  as  sea. 


2. 


Transforms  and  Measured  Radiances 


2 . 1  Production  of  the  Transforms 

A  general  capability  has  been  developed  which  relates  any  radiance, 
for  which  the  transmission  function  is  specified,  to  parameterizations  of 
the  atmospheric  thermal  structure.  This  program  represents  a  merging  of 
items  (1)  and  (2)  listed  under  Section  ]  .1  above,  but  with  added 
generalizations  and  extensions. 

The  program  permits  arbitrary  specification  of  (a)  the  standard  levels 
which  define  the  atmospheric  mass,  and  thermal  structure,  parameterization, 
(b)  a  standard  temperature  structure  in  terms  of  values  for  the  degrees  of 
freedom  of  the  modelled  thermal  structure,  and  (c)  the  transmission  function, 
for  any  channel  and  nadir  angle,  tabled  in  terms  of  a  specified  increment 
in  p*. 

The  method  of  integrating  the  transmission  function  for  the  standard 
thermal  structure.  the  linearization  of  the  Planck  function  about  the  standard 
thermal  structure,  and  the  formation  of  the  transforms  are  detailed  in  the 
Task  Two  Report  of  Phase  One. 

For  development  purposes  we  chose  an  18-level,  16-layer  mass- 
structure  model.  The  specified  standard  pressure  levels  are  listed  in  the 
first  column  of  Table  1 . 

The  arbitrary  standard  temperature  structure  that  we  have  selected 
is  also  specified  in  Table  1.  The  integration  of  the  radiative  transfer 
equation  is  performed  in  terms  of  the  layers  in  which  temperature  is  modelled 
to  be  linear  in  p*.  The  top,  layer  interface,  and  bottom  temperatures  of  the 
assumed  standard  thermal  structure  are  listed  in  the  last  column  of  Table  1. 
This  standard  atmosphere  is  plotted  on  a  Sigmagram  in  Fig.  1.  The  ground, 
or  sea  surface,  .emperature  is  specified  to  be  the  same  as  that  for  1000  mb. 


Table  1 


Specification  of  the  arbitrary  standard  atmosphere  may  be 
expressed  in  any  of  the  following  terms  which  are  mutually 
consistent  with  the  18-level,  1  6-laycr  mass-structure  model 


By  Heights 


Model-  Specified 

Standard  Arbitrary 

Pressure  Standard 

Levels:  Heights: 

(mb)  (meters) 


By  Thickness 
and 
Layer 
Stabilities 


Layer:  Stability: 

(mb)  (102m2sec“2) 


By  Temperature 
at  Top 

Interface  and 
Base 

Level  Temperature 
(mb)  (c.  K) 


0 

79,000 

(Z300  "  Z1000  ~ 

9 , 063  meters) 

0 

249.8 

4 

37,100 

0-5.5 

20.49 

5.5 

224.9 

7 

33,400 

5.5  -  8.5 

16.40 

8.5 

228.2 

10 

31,024 

8.5  -  15 

20.04 

15 

225.3 

20 

26,451 

15  -  25 

21.08 

25 

220.3 

30 

23,831 

25  -  40 

19.61 

40 

217.6 

50 

20,574 

40  -  60 

19.62 

60 

215.0 

70 

18,455 

60  -  85 

19.11 

85 

213.1 

100 

16,226 

85  -  125 

16.79 

125 

214.1 

150 

13,681 

125  -  175 

14.11 

175 

218.3 

200 

11,843 

175  -  225 

13.57 

225 

222.3 

250 

10,391 

225  -  275 

12.14 

275 

226.6 

300 

9,180 

275  -  350 

7.27 

350 

236.5 

400 

7,191 

350  -  450 

6.29 

450 

248.4 

500 

5,571 

450  -  COO 

7.28 

GOO 

2  62 . 0 

700 

2,998 

600  -  775 

8.72 

775 

273.9 

850 

1,445 

775-1000 

13.95 

1000 

281  .7 

1000  117 


T  (ground)  =T  (1000  mb) 
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The  definitions  and  relationships  may  be  identified  as  follows: 

E  (v  ,6)  denotes  the  radiance  measured  by  channel  v  at  nadir  angle  B  . 

Eg  (v  ,0)  denotes  the  standard  radiance  composed  of  a  ground,  or  sea- 
surface,  contribution  and  an  atmospheric  contribution.  The 
program  produces  standard  values  for  both  contributions. 

The  radiance  anomaly  is  defined  by 

e  (V  ,B)  f  e(i/,6)  -  Es(v,6) 

The  radiance  anomaly  is  related  to  the  temperature  anomalies,  T,  at  the  ground 
and  at  the  top,  layer  interfaces,  and  bottom  of  the  modelled  thermal 
structure: 

Mt'.S)  =  Cn(v,  B)3  +  £  C  (i/,0)3:  .  (1) 

1  n  n 

n 

The  program  calculates  the  values  of  these  coefficients,  Cq  for  the  ground, 
and  the  Cn's  for  the  column  levels.  These  column  levels  are  identified  in 
Table  2  under  the  heading  C  (I)  . 

By  virtue  of  the  intentional  linear  properties  of  the  mass-structure 
model  the  radiance  anomaly  may  also  be  expressed  in  alternate  terms  of 
the  thermal-structure  anomaly: 

A  linear  combination  of  the  height  anomalies  of  the  standard 
pressure  levels  as  identified  in  Table  2  under  the  heading  CA(I). 

The  corresponding  combination  coefficients  always  sum  to  zero 
because  the  thermal  structure  is  not  related  to  the  absolute  of 
the  height  structure. 
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Table  2 


the  columns  or  coefficients,  listed  for  each  frequency  akd  ancle, 

ARE  LINEAR  COMBINATION  UOEI F l C I E NTS  FqR  EXPRESSING  THE  RADIANCE 


ANOMALY 

IN  TERMS  OF 

the  following 

SETS  OF  PARAMETER 

anomalies? 

300^1000  MB 

I 

cm 

cam  ) 

CO  f  I  ) 

CCf  I  ) 

1 

T  (  0  ) 

Z<  0 ) 

Z(  0)«z(  4) 

SIGKA16 

2 

T  (  5) 

Z<  A  ) 

Z(  4  )*Z<  7) 

SJGMA15 

3 

T<  8) 

Zt  7) 

Z(  7)*z(  10) 

SIGMA14 

4 

T  (  15) 

Zv  10) 

Z  (  10)*Z<  *0) 

SIGM13 

5 

T  (  25) 

Z(  20) 

Z(  2Q)*Z<  <30) 

S1GWA12 

6 

T  (  40) 

7A  30) 

Z(  30  )  *Z  (  !>0) 

SIGKA11 

7 

T  (  60) 

Zl  50) 

Z<  50  )  »Z I  7o) 

SIGWA10 

8 

T  (  85) 

Z<  70) 

Z(  70 ) «Z ( 10 0 ) 

SIGMA  9 

9 

T ( 125 ) 

ZU00  ) 

Z(1Q0  )«Z(1S0 > 

SIGMA  a 

10 

T  < 1 75 ) 

Z ( 150 ) 

Z 1 150 ) *Z (200  > 

SIGMA  7 

11 

T (225 ) 

ZI2Q0) 

Z(200)*Z(290) 

SIGMA  6 

12 

T ( 275 ) 

Z  C  250  ) 

Z ( 250  ) *Z ( 30 0  ) 

SIGMA  5 

13 

T I 350  ) 

Z  l  3  C  0  ) 

Z(3CO)-*Z<400  ) 

sigma  4 

14 

T ( 450  ) 

Z  (  4  o  o  ) 

Z(400)*ZI5U0) 

SIGMA  3 

15 

T ( 60 0  ) 

Z  t  5g  0  ) 

Z(500  )*Z(/00  ) 

SIGMA  2 

16 

T<775> 

Z  I  7  o  0  ) 

Z(7oo)-z(as0) 

SIGMA  1 

17 

T  <  1 D  0  0  ) 

Z<«50) 

Z(850)*ZC1000) 

W 

18 

Z(1000) 

Z(iooo) 
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A  linear  combination  of  the  thickness  anomalies  of  layers,  defined 
by  the  succession  of  standard  pressure  levels,  as  identified  in 
Table  2  under  the  heading  CB(I). 

A  linear  combination  of  the  anomalies  of  a  single  thickness,  H  ,  and 
of  the  static  stabilities  of  each  of  the  16  layers,  as  identified  in 
Table  2  under  the  heading  CC(I).  We  have  chosen  the  thickness 
from  1000  to  300  mb;  most  of  the  variance  of  the  atmospheric  thermal 
structure  is  resolved  by  specification  of  this  thickness. 

We  consider  these  to  be  the  major  expressions.  Others  can  be  designed  within 
the  linearities  of  the  mass-structure  model. 


We  were  provided  with  the  transmission  functions  for  an  operational 
Air  Force  VTPR  sounder,  for  a  set  of  channels  at  a  set  of  nadir  angles.  Each 
transmission  function  was  tabled  in  69  equal  increments  in  p*,  from  1000  mb 
to  0.01  mb.  We  produced  a  complete  set  of  standard  radiances  and  anomaly 
combination  coefficients  for  this  sounder.  A  sampling  is  here  included  in 
Tables  3,  A,  B  and  C.  The  legend  corresponds  to  Table  2. 
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Table  3  A 

Transformation  Table  for  Channel  608.4  angles  25°  and  30° 

FREQUENCY  =  6AH.40  ANGwfc  a  25,00 

ASSOCIATED  PRESSURE l  P=  11,96  P**KAPPa=  2,031933 

STANDARD  RADIANCE  s  49.99293//1P  =  0,000000000  ♦  49,992937710 

radiance  anomaly j  gce  o.ououuqpoo  x  t ( r; > 


3QU-1000  M=» 


I 

C  (  I  ) 

CAl  I  ) 

C3(I) 

ecu ) 

1 

,122530659 

.001188342 

,001188542 

,001494063 

2 

,274365937 

.011421/84 

,012610326 

,000391261 

3 

,091101754 

-  ,008396365 

,004213960 

',00061)8934 

4 

,112219104 

,002955201 

,007169161 

,000655156 

5 

,111508215 

,000807216 

,0079/6380 

,000694772 

6 

,096231677 

-.000352808 

,  0  P  7  623577 

,000657161 

7 

,070346954 

-.00119997/ 

,006423600 

’,000588068 

8 

,049627105 

-.001471/89 

,004951811 

‘,000649513 

9 

,025381471 

-  .002599823 

,002351906 

,000539430 

10 

,006691675 

~,  001799339 

,000552646 

,000376215 

11 

,001112478 

-,000441603 

,000111041 

,000282352 

12 

,000128754 

-,000101162 

,000009879 

,00030/824 

13 

0,000000000 

-,000010645 

-.00001)0766 

,00024697/ 

14 

0,000000000 

,000000922 

,000000155 

,000192570 

15 

0,000000000 

-,00000016b 

-,000000013 

7000096329 

16 

0,000000000 

,00000001/ 

.000000QU4 

',  000030902 

17 

0,000000000 

-  .OOOOOOU04 

-  ,  000000001 

,009102944 

18 

.OOOOUOUOl 

,000000000 

FREQUENCY  =  668,40  ANGLE  *  30,00 

ASSOCIATED  PHESSURE  J  P  =  11,55  P**KAPPA=  2,011899 

STANDARD  RADJANCE  =  50.140/62/27  =  0,000000000  *  50,140/6272? 

RADIANCE  ANOMALYJ  GCb  0,000000000  X  T(G> 


30U-1000  M8 


I 

cu  > 

CAU) 

C8U  > 

ecu ) 

1 

,126234032 

.001225849 

,001225849 

,001539225 

2 

,277934210 

,091478368 

,01143332/ 

,012659176 

,000396349 

3 

-.008429966 

,004229210 

.000615378 

4 

,112689724 

,002975152 

,007204362 

,000660901 

5 

,111295755 

,  000759946 

,  Q 0796430a 

',  000699041 

6 

,095373879 

-,000407476 

,007556832 

‘,000659272 

7 

,069118783 

-,001243223 

,006313609 

,000588278 

L 

,047906514 

-.001528672 

,004784937 

,000647533 

9 

,023797600 

-.002589390 

,002195547 

,000536102 

10 

,006117430 

-,00172845b 

,000467091 

.U00373276 

il 

,000947531 

-.000371371 

,000095720 

,000280040 

12 

,000099008 

-.000090122 

,000005598 

,000305281 

13 

0,000000000 

-.000006032 

•,000000434 

,’U  0024493/ 

14 

0,000000000 

.000000522 

,0O00UOQ88 

V  0  0  0190  979 

15 

0  .OUQOOOOOO 

-.000000095 

-  ,  000000007 

,000095534 

16 

0  ,  00000001)0 

.oooooooiu 

,000000002 

,000030646 

17 

0,000000000 

-,  000000003 

■'OOQOQPQOO 

V 009027743 

18 

,000000000 

,000000000 
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Table  3  B 

Transformation  Table  for  Channel  677  angles  55°  and  57.5° 

FREQUENCY  *  677,00  ANGUfc  «  55,00 

ASSOCIATED  PRESSURE  I  Pa  22,6©  P**KAPPA=  8,439937 
STANDARD  RADIANCE  a  46.07873/826  =  0 ,000000000  *  46 , 078  737.8  2* 

RADIANCE  ANOHALYI  GCx  O.OUOOOQOOO  X  T ( G > 


300-1000  MB 


I 

C  ( 1 ) 

CAt  I  > 

cBtn 

CCU) 

1 

,054038693 

,0005l8b6» 

,000518669 

;U006539l7 

2 

,142336372 

.006363413 

,006882Q8l 

,000202979 

3 

,092283641 

-,00264424b 

,004237336 

,000379874 

4 

,120699380 

,003326148 

,007563984 

1000492649 

9 

,140820787 

,002472131 

,010036115 

.UQ0609732 

6 

,138937842 

,0009/9630 

,011015745 

,000649398 

7 

,111795766 

•,00068058a 

,010365157 

7000629955 

8 

,082335241 

-,002062612 

,008302943 

7000732891 

9 

,040031277 

-.004606014 

,003696531 

7000621408 

10 

,008626259 

-.003165470 

,000531061 

700Q43393V 

11 

,001017929 

-.000436153 

,000091909 

7U00323517 

12 

,000156742 

-.000075261 

,000019547 

,000354839 

13 

,000003066 

000020 63  V 

-  .  000000393 

,000284742 

14 

0,000000000 

,000001194 

,000000301 

.00022201? 

15 

0,000000000 

-.000000218 

.,000000317 

,000111059 

16 

0,000000000 
o', ooooooooo 

,000000022 

.000000005 

,00003562/ 

17 

-,000000006 

..OOOOOQOOl 

.010494845 

IB 

,000000001 

.uoonooooo 

FREQUENCY  ■  677,00  ANGLE  a  57,50 

ASSOCIATED  PRESSUREl  P»  2l,7{)  P**KA9PA=  2.409060 

STANDARD  RADIANCE  »  46,265691395  «  O.OOQOQOOOO  ♦  46,265691395 


RADIANCE  ANOHALYI  GCa  0,000000003  X  T(G> 


I 

CU> 

CAM  ) 

cam 

300-1000  MB 
CCU  ) 

1 

,056800186 

,000545745 

,000345745 

,00069258V 

2 

,147597779 

,006550803 

,007096540 

1000210432 

3 

,084451583 

-,002734817 

,004361731 

,000392320 

4 

,123156720 

,003363943 

.007723074 

,000506311 

5 

,142251483 

,002429075 

,010134749 

;U0062284V 

6 

,138629315 

,000842882 

,010997631 

,’000658970 

7 

,109726187 

-,000320112 

,  0101 7 7 319 

00C6349-46 

8 

,070462247 

-.00224547J. 

,007932047 

V00073H941 

9 

,036091009 

• ,  0046371*18 

,003294229 

,000617044 

10 

,007102529 

-.002911696 

,000332533 

,000429686 

11 

,000756239 

-,00030983V 

,000072094 

000322133 

12 

,000127179 

-,000055952 

,000016741 

,000351178 

13 

,000003066 

-.00001750^ 

.,000000767 

JO0Q281764 

14 

0,000000000 

,000000923 

,000000159 

,000219693 

15 

0,000000000 

-.000000168 

.,000000013 

,00010989/ 

16 

0,000000000 

,000000017 

,000000904 

,000035254 

17 

0,000000000 

-,000000004 

.,000000901 

,U 1038503 7 

IB 

.oooooouox 
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Table  3  C 

Transformation  Table  for  Channel  746.5  angles  0°  and  10° 

FREQUENCY  *  746,50  ANGLE  >  0,00 

ASSOCIATED  PRESSURE  I  P«  700,4*  P**KAPPAp  *,500743 
STANDARD  RADIANCE  =  09,219097408  ■  45,277022937  *  43,94207447 t 
RADIANCE  ANOMALYJ  GC.  ,62659*317  X  T(Q> 


300*1000  MO 


I 

CM) 

CAM) 

CB(I) 

CC<  I  ) 

1 

,001848667 

,000017840 

,000017840 

,000022542 

2 

,004523313 

•  0  00194249 

.000212388 

.U0000645U 

3 

,002490520 

*,000089051 

,000123338 

,000011843 

4 

,003827063 

,000113334 

,000236672 

,000015466 

5 

,005249576 

,000114893 

,000351567 

,000020412 

6 

,006834653 

,000176347 

,000527914 

'.000024652 

7 

,000232029 

,000211913 

,000739828 

,00002*520 

6 

,011723397 

.000281151 

,001020979 

7000043353 

9 

,016459664 

,000515452 

,001536*31 

V000052508 

10 

,019005922 

,000707007 

,002243*30 

,003179910 

;oooo>2iii 

u 

,020837654 

,000936471 

’,000053163 

12 

,032492660 

.001187951 

,004367860 

7000062636 

13 

,058678837 

,003492/37 

,007860597 

-.000003149 

14 

,103728607 

,004295510 

,012156107 

., 000054637 

15 

,153671655 

,006944593 

,019100700 

-7000048590 

16 

,190206422 

,009419203 

,028519983 

-.000009330 

17 

,083367789 

*,006774252 

,021745730 

,019973231 

18 

-.021745/3U 

7000000000 

FREQUENCY  »  746,50  ANGLE  s  10,00 

ASSOCIATED  PRESSURE:  P=  690,19  P**KAPPA=  *,494772 
STANDARD  RADIANCE  *  89,046456921  ■  44,9077343*7  *  44,130732604 
RADIANCE  ANOMALY!  GC*  ,621483606  X  T(G> 


300*1000  MS 


l 

CM  ) 

CAM  ) 

C9U) 

CC<  I ) 

1 

,001850287 

,000017828 

,00001782b 

.000022561 

2 

,004612321 

,000201269 

,000219097 

',000006577 

3 

,002490520 

- , 0  0  00  99265 

,000119032 

,000011988 

4 

,003837009 

.000120715 

,000240540 

7009015676 

5 

,005340843 

, 000117303 

,001)357)31 

7000020720 

6 

,006940646 

,000178840 

,00  J 536691 

,000025028 

7 

,008306311 

,000200970 

,  000  ’45661 

7000028896 

8 

,011032381 

,000204218 

,001029*79 

7000043066 

9 

,016639976 

,000523728 

,001553606 

7000053113 

10 

,019197358 

,00071372V 

, 002267*35 

,000052688 

11 

,020979170 

.000931/78 

,003199113 

7000053691 

12 

,032734340 

,001202791 

,004401905 

7000063325 

13 

,059095734 

,003515105 

,007917910 

•7000002621 

14 

,104337174 

,00431978* 

,012236793 

•7000054148 

15 

,154302546 

,006942040 

,019178*33 

.,000048235 

16 

,190744654 

,009433719 

,026612552 

•7000009Q94 

17 

,0*3429948 

-,006861197 

,  021751*55 

7019959012 

18 

-.021751355 

-20- 

7000000000 

% 

2.2  Normalization  of  the  Radiance  Anomalies 


The  purpose  of  the  present  section  is  to  develop  a  feeling  for  the 
resolution  powers  of  this  Air  Torce  sounder.  The  combination  coefficients 
0(1),  which  apply  to  the  temperature  anomalies,  are  all  positive — suggesting 
the  concept  of  a  smoothing  operation.  A  radiance  measure  is,  in  fact,  a 
smoother;  it  has  such  a  broad-in-depth  energy-source  function. 

Application  of  the  concept  of  a  smoother  requires  normalization  of 
the  combination  coefficients.  Unless  already  normalized, 

£cn  t  1 

Normalization,  i.e.,  scaling  the  radiance  anomaly  as  a  smoothed  temperature 
anomaly,  is  accomplished  by  dividing  the  radiance  measure  by  £  Cn  .  This 
requirement  is  obvious  when  applying  Eq.  (1)  to  a  temperature  profile  of 
constant  anomaly,  T n  =  T  .  we  obtain 

«/Ecn  =  3 

If  we  are  to  consider  a  normalized  radiance  as  a  value  in  a  smoothed 
temperature  profile  then  we  must  also  be  able  to  associate  the  value  with  some 
level  in  the  profile.  In  the  preceding  we  have  explained  normalization  on  the 
basis  of  a  constant  temperature-anomaly  profile.  We  now  use,  as  basis  for 
the  definition  of  an  associated  pressure  level,  a  temperature  anomaly  profile 
which  is  linear  in  p*,  including  the  surface  value.  Such  a  profile  is  not 
altered  by  a  smoothing  operator  which  is  symmetric  in  pK,  when  applied 
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at  any  level,  p  .  On  this  basis  wo  define  the  associated  pressure  level  as 
that  level  where 

=  rr 
c 

n 

for  If  linear  in  p*: 


rr  -  /3(afp«)  .  (3) 

The  selection  of  interfaces,  pn  ,  in  the  modelling,  is  immaterial  for 
this  determination.  Any  segment  of  Eq.  (3)  is  linear  in  p*  and  the  segments 
join  up  at  the  interfaces  to  form  a  complete  linear  profile. 

Substitution  of  Eq .  (3)  into  Eq.  (2)  yields  the  associated  pressure 

level 


P 


Ec  p* 

n  n 


1 A 


(4) 


Note  that  Eq.  (4)  is  independent  of  the  a  and  of  Eq.  (3).  The 
implication  is  not  that  all  radiances  are  equivalent  to  symmetric  smoothing 
operators .  An  asymmetric  smoother  applied  at  any  level  in  a  linear  profile 
produces  an  increase  or  decrease  at  that  level.  Our  definition  of  the 
associated  pressure  level  accordingly  involves  an  upward  or  downward 
displacement  of  the  level  of  application  of  an  equivalent  smoothing 
operator.  This  shift  is  pronounced  for  radiances  having  large  surface 
contributions . 
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Calculation  of  this  associated  pressure  level,  for  each  channel 
and  nadir  angle,  has  been  added  to  the  general  program  which  produces  the 
combination  coefficients.  These  values  are  found  on  the  second  line  of 
Tables  3  A,  B  and  C,  expressed  both  in  mb  and  in  units  of  p*.  The 
normalizing  factor  is  obtained  by  summing  column  C(I)  and  adding  the 
coefficient  for  the  ground  anomaly  contribution  found  on  the  fourth  line  of 
those  Tables . 

Table  4  summarizes  the  standard  radiance,  the  normalizing  factor, 
and  the  associated  pressure  level  in  mb  and  in  units  of  p*,  for  each  channel 
and  nadir  angle  of  the  Air  Force  VTPR  sounder. 
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Table  4 

Summary  of  Channels  for  Specified  Zenith  Angles 


FRED 

ANU 

97450, NAD, 

N|RH, FACTOR 

P 

(••KAI 

834,0 

0,00 

V9, 036490011 

1,51840*021 
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0,00 

40,600947040 
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2(075, 

668,4 

10,00 

49. 728/60202 

,9906*7*67 

19,74 

2 1 0  69, 

688,4 

20,00 

49,6/6002207 

,900225510 

12,99 

2 1  0  4  7' 

668,4 

29,00 

49,942437710 

,901445712 

11.90 

2l  0  31' 
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30,00 

50.140/6J727 
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11.91 
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38,00 

50,3234241X7 

,9449141*9 

11.07 

1,907, 
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40,00 

50,847140630 
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10.91 

1(988' 
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42,90 

90,0/6/64095 

,900414*02 

11.21 

1(»42 

668,4 

49.00 

50,010400704 

,9701)0*45 

0.00 

ll»24 

668,4 

47,90 

80,9/5/09137 

,971737685 

9.90 

1(904 

668,4 

90,00 

51.140402919 

,9719398*4 

0.10 

1,001 

668,4 

92.90 

91,339265709 

,9755201/9 

0.79 

1(961 

668,4 

99,00 

51,549209174 

,977/10215 

0.29 

1,916 

668,4 
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51,703010239 

,960130719 

7,91 

1,009 

668,4 

60,00 

52,042167700 

,902022689 

7,51 

1(700 

677,0 

0.00 

44,7l6»6j«02 

,907950930 

33,00 

2.U6 

677,0 

10,00 

44,747762727 

,900290355 

32,70 

2,709 

677,0 

20,00 
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31.00 
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677,0 

29,00 

44,923/43096 
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30.9) 

2(069 

677, ( 

30,00 

49,024500312 
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21,99 
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677,0 

38,00 

49,152/00011 
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29,00 
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677, t 

40,00 

49,314175292 
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27,02 

2(580. 

677,0 

42,90 

49,400551241 

,9154*3077 

20/91 

2,561 
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49,00 

49,514193909 
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29.19 
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47,80 
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21,39 
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677,0 
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24,92 

2,494 

677,0 

92.80 
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01,71 
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79,97 

3/446 
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,896102813 

0),89 
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0*.  1* 
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0,00 

59,641*29610 
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2>1.19 
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707,8 

10,00 
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20,03 
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30,00 
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40,00 
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0,624 
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45,00 
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7*7,5 

97,90 
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1,114110/90 

333,27 
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6,900 

74...  8 

10,00 

89.044456921 

1,34020*903 

099,10 

6(494 

746,8 

20,00 

08,511503439 

1,34337032* 

*»1.12 

0(479 
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2 . 3  Radiance  Anomaly  Expectations 


Based  on  the  transformation  tables,  we  calculated  the  normalized 
radiance  anomalies  to  be  expected  for  a  sampling  of  specified  temperature- 
anomaly  profiles.  The  purpose  of  these  samples  is  to  evaluate  our 
contention  that  the  normalized  radiance  anomalies,  each  at  associated 
pressure  levels,  can,  to  a  degree,  be  regarded  as  equivalent  smoothed 
tcmperature-anomaly-profile  values . 

The  profile  is  specified  in  terms  of  the  temperature  anomaly  of  the 
ground  and  at  each  of  the  interfaces  identified  in  column  one,  C(f),  of 
Table  2.  The  expected  normalized  radiance  anomaly  was  calculated  by 
linearly  combining  the  specified  temperature  anomalies  v/ith  the  coefficient 
sets  corresponding  to  each  channel  and  nadir  angle,  and  then  dividing  the 
summation  by  the  sum  of  these  coefficients  in  order  to  normalize  to 
temperature  scale.  For  each  specified  profile  a  plot  was  prepared  which 
shows  the  specified  temperature-anomaly  profile  by  connected  linear- 
in-pK  segments  in  accordance  with  the  mass-structure  model;  the  calculated 
normalized  radiance  anomaly,  for  each  channel  and  each  nadir  angle,  is 
shown  by  a  dot  plotted  at  associated  pressure  level. 

Samples  are  given  as  Figs.  2,  3,  4,  5,  and  6.  The  vertical  scale 
is  linear  in  p*,  rescaled  in  mb  on  the  right.  The  horizontal  scale  is  linear 
in  temperature  anomaly;  the  calibration  is  irrelevant.  The  profile  connects 
horizontally,  at  1000  mb,  to  the  ground  temperature  anomaly  if  it  differs 
from  the  air  temperature  anomaly  (See,  e.g.,  Fig.  3). 

Figure  2  is  a  demonstration  showing  the  general  character  of  the 
sounder.  Each  of  the  six  strings  of  dots  corresponds  to  one  of  the  six 
channels  as  identified  in  Table  4. 

Figure  3  demonstrates  the  ground  contribution. 


-25- 


Figure  4,  depicting  a  cold  anomaly  ranging  from  450  mb  to  85  mb. 


demonstrates  that  the  normalized  radiance  anomalies  cannot  be  strictly 
interpreted  as  smoothed  temperature-anomaly-profile  values .  Note  the 
discontinuities  between  channels.  These  discontinuities  do  not  relate 
to  channel  calibrations,  but  rather  to  the  differential  character  of  channel 
source  functions.  Figure  5  shows  that  the  shift  between  channels  can 
also  be  of  the  other  sign.  This  disparity  in  channel-to-channel  smoothness 
is  associated  with  enhanced  vertical  discrimination  powers  of  the  sounder. 
By  the  same  token,  the  angle-to-angle  smoothness  for  any  one  channel 
demonstrates  that  little  is  added  as  to  vertical  discrimination  powers  by 
views  over  a  range  of  zenith  angles. 

Figure  6  demonstrates  the  concept  of  the  normalized  radiance 
anomaly  and  associated  pressure  level.  This  anomaly  profile  is  linear 
in  pK.  All  dots  fall  on  the  profile  by  definition. 

The  surface  contribution  can  be  removed  from  each  channel  radiance 
by  sacrificing  one  channel  for  that  purpose.  For  example,  the  window 
channel  may  be  used  for  that  purpose.  However,  in  our  calculated  profile 
examples  this  is  accomplished  by  merely  setting  the  coefficient  for  the 
ground  contribution  to  zero.  The  normalizing  factor,  and  associated 
pressure  for  each  radiance,  change  according  to  definition.  Figures  7,  8 
and  9  are  a  repeat  of  profiles  shown  in  Figs.  2,  4  and  5  respectively, 
without  ground  contribution  in  the  calculated  radiances. 

These  examples,  with  and  without  ground  contribution,  demonstrate 
the  following:  the  angle-to-angle  normalized-radiance-anomaly  continuity 
for  any  one  channel  is  pronounced  and  can  be  interpreted  as  a  smoothing 
of  the  temperature-anomaly  profile  at  associated  pressure  levels.  Whereas 
a  coarse  continuity  in  normalized  radiance  anomaly  exists  from  one  channel 
to  the  next,  strict  continuity  is  not  to  be  imposed  in  processing  true 
radiances,  or  in  relative  calibrations. 
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Fig.  3  A  specified  temperature  anomaly  profile  and  plotted  values  of 
calculated  normalized  radiance  anomalies. 
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Fig.  5  A  specified  temperature  anomaly  profile  and  plotted  values  of 
calculated  normalized  radiance  anomalies. 
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Fig.  7  A  specified  temperature  anomaly  profile  and  plotted  values  of 
calculated  normalized  radiance  anomalies  without  ground 
contribution. 


Fig.  8  A  specified  temperature  anomaly  profile  and  plotted  values  of 
calculated  normalized  radiance  anomalies  without  ground 
contribution. 
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Fig.  9  A  specified  temperature  anomaly  profile  and  plotted  values  of 
calculated  normalized  radiance  anomalies  without  ground 
contribution. 
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2.4  Processing  of  Real  Satellite  Data 


We  were  provided  with  a  sample  of  raw  radiances  measured  by  the 
Air  Force  VTPR  sounder  in  orbits  of  March  26,  1974.  These  data  include  full 
side-to-side  scans  with  the  position  of  each  spot  given  in  latitude  and 
longitude,  and  with  nadir  angle,  and  time  data. 

The  specified  nadir  angles  are  not  limited  to  the  discrete  values 
for  which  we  have  generated  tables.  For  nadir  angles  which  fall  between 
our  tables,  linear  interpolation  is  used  for  tabled  quantities. 

We  have  processed  each  complete  side-to-side  scan  into  a  cross 
section.  Standard  radiances  are  removed  to  form  radiance  anomalies  which 
are  then  normalized  for  plotting  at  associated  pressure  levels.  The  decimal 
point  marks  the  cross-section  position  of  each  normalized  radiance 
anomaly.  The  vertical  coordinate  of  the  cross  section  is  linear  in  p*. 
Sample  cross  sections  begin  with  Fig.  10. 

The  horizontal  coordinate  of  the  cross  section  is  calibrated  in  units 
of  standard  grid  length  measured  along  the  scan  base  line.  For 
each  spot  viewed,  the  grid  locations  are  given  below  the  base  line,  with 
I  above  J  .  These  grid  coordinates  are  for  the  standard  63x63  grid  array 
for  the  northern-hemisphere  polar-stereographic  projection.* 

The  purpose  of  these  cross  sections  is  a  graphical  rendering  to 
facilitate  subjective  examination  of  the  character  of  the  radiances  and 
the  nature  of  cloud  contamination.  In  a  cloud-free  section,  over  open 
ocean,  the  normalized-radiance-anomaly  values  for  any  one  channel, 
strung  across  the  section,  can  generally  be  expected  to  be  relatively 
smooth  (See,  e.g.,  Fig.  10).  Any  lack  of  smoothness  generally  implies 
cloud  discontinuities,  or  abrupt  changes  in  surface  temperature.  In  both 
cases,  cloud  or  ground  induced,  the  discontinuities  become  less  pronounced 
for  channels  with  higher  (i . e . ,  smaller  pressure)  associated  levels, 

* 

As  shown  in  Fig.  13,  page  59. 
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reflecting  a  decrease  in  the  contribution  of  lower  levels  to  these  channels . 
These  characteristics  are  demonstrated  in  Fig.  12. 


About  eighty  of  these  cross  sections  were  produced.  On  the  whole, 
the  radiance  values  appear  credible.  The  discontinuities  from  channel  to 
channel  generally  appear  to  be  small  and  acceptable.  In  some  cases — such 
as  in  Figs.  10  and  11 --the  discontinuities  between  lower  channels  appear 
to  be  rather  prom  unced.  However  we  are  not  as  yet  in  a  position  to 
comment  on  relative  and/or  absolute  calibration  of  channels. 

We  have  studied  these  cross  sections  and  other  graphical 
formats  in  order  to  evaluate  cloud-top  models  and  in  order  to  develop  our 
approach  for  diagnosing  the  radiances  for  their  clear-column  components. 


Fig.  10  Cross 


N 

f'J 


Fig.  12  Cross  Section  with  Ground  Contributions  and  Window  Channel 


3. 


Development  of  CLRX 


3 . 1  Introduction 

The  major  development  under  Phase  Two  is  a  basic  capability  for 
diagnosing  the  clear-column  radiance  components  from  cloud -contaminated 
VTPR  radiances.  The  basic  scheme,  designed  for  applicability  to  the  Air 
Force  VTPR  system,  processes  one  complete  scan*  at  a  time.  We  refer  to 
this  programmed  capability  as  Program  CLRX. 

The  new  procedure  is  based  on  assessing  various  forms  of  information 
inherent  in  the  radiances,  and  blending  all  this  information,  simultaneously 
for  the  full  scan,  to  produce  the  optimum  estimates  of  clear-colum  radiances 
for  each  spot.  All  input  information  and  modelling  statements  are  weighted 
in  inverse  proportion  to  the  associated  estimated  variances.  The  produced 
optimum  estimates  of  the  clear-column  radiances  are  those  which  minimize 
the  total  disparity  in  accommodating  all  weighted  elements  of  information. 
The  procedure  also  produces  estimates  c i  the  associated  reliability  for  each 
produced  clear-column  radiance  element;  i.e. ,  for  each  channel  and  scan 
spot.  This  new  procedure  is  another  application  of  our  Fields  by  Information 
Blending  (FIB)  methodology,  introducing  additional  concepts  and 
formulations . 

The  basic  version  exploits  the  following  elements  of  information: 

(a)  A  tentative  non-zero  reliability**  is  assigned  to  each  VTPR  radiance 
element  which  is  estimated  to  be  close  to  the  clear-column  component 


* 

The  VTPR  scans  of  NOAA  spacecraft  are  more  closely  spaced  along  the 
orbital  path,  affording  useful  information  redundancies.  This  redundancy 
warrants  exploitation  by  analyzing  differences  between  measured  spots 
not  just  within  one  scan  but  also  between  adjoining  scans. 

★  V 

Ir  the  FIB  methodology  reliability  is  generally  defined  as  the  inverse  of 
the  associated  error  variance . 
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for  that  element.  This  assessment  is  based  on  the  variance  character 
for  spots  and  channels  in  the  ambience  of  the  element. 

(b)  An  information  statement  which  states  that  the  clear-column  component 
for  each  channel  does  not  vary  abruptly  from  spot  to  spot  across  the 
scan.  This  smoothness  condition  can  generally  be  imposed  on  all 
channels  in  oceanic  regions.  The  basic  version  is  designed  for 
application  in  oceanic  regions. 

(c)  A  model  for  interpreting  cloud  contamination.  The  interpretation  model 
is  based  on  the  following  consideration:  For  any  pair  of  spots,  in 
proximity  to  each  other,  the  largest  difference  between  them — of 
significance  to  the  radiation — is  generally  the  amount  of  cloud  of 

a  single  cloud-top  level.  Even  a  cursory  inspection  of  satellite 
videographs  confirms  this  impression.  Based  on  this  interpretation, 
large  differences  in  radiance  between  nearby  spots  produce  estimates 
of  the  correction  axis  (i.e. ,  unit  vectors)  for  adjusting  the  measured 
radiances  toward  their  clear-column  components. 

(d)  Independent  Sea-Surface  Temperature  (SST)  information.  This  new 
capability  does  not  involve  any  categorical  use  of  independent  SST 
information.  However,  any  available  estimates  of  SST,  and/or 

SST  gradient  along  a  scan,  can  be  exploited  to  enhance  the  resolution 
of  the  clear  radiance  components.  Other  independent  information, 
especially  radiance  estimates  derived  from  preliminary  upper-air 
analyses  and  extrapolation,  can  be  exploited  by  this  capability. 
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The  application  of  the  FIB  methodology  includes  the  following 
component  operations: 

•  The  input  of  VTPR  radiances,  and  the  derivation  and  weighting  of  all 
elements  of  information,  is  called  the  assembly  operation.  Assembly 
includes  input  of  independent  sources  of  information  such  as  transformed 
sea-surface  temperature  estimates  and  difference  estimates  along  the 
scan  base. 

•  Assembly  is  followed  by  the  blending  operation  which  may  be  performed 
by  iterative  solution  or  by  explicit  matrix  inversion  schemes.  The 
blending  operation  produces  the  resultant  of  all  assembled  information: 
the  clear-column  radiance  elements  and  their  associated  reliabilities. 

•  The  reevaluation  and  error-checking  operation  follows  the  blending. 

All  tentatively  assigned  weights  are  subjected  to  reevaluation  in 
comparison  to  the  corresponding  implications  of  all  information. 

•  Recycling:  The  blending  is  repeated  using  reevaluated  information 
statements.  This  is  followed  by  a  second  reevaluation  and  a  final 
blending. 

The  resultant  clear-column  estimates  will  generally  be  of  uneven 
accuracy  as  indicated  by  the  distribution  of  associated  reliability  estimates. 
There  will  be  great  variations  in  associated  reliabilities  within  scans, 
varying  from  channel  to  channel  and  spot  to  spot,  and  there  will  be  large 
variations  from  scan  to  scan.  The  associated  distribution  of  reliability 
estimates,  produced  for  each  scan,  will  depend  on  the  quality  of  the  input 
information  and  on  the  character  of  the  scan:  the  distribution  of  cloudiness 
and  the  degree  of  spot-to-spot  compliance  with  the  imposed  interpretation 
model. 

The  associated  reliabilities  can  be  used  for  categorical  acceptance 
or  rejection  of  the  individual  diagnosed  clear-column  radiance  components. 
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Properly  tuned,  the  associated  reliabilities  can  be  more  fully  exploited  as 
weights  for  the  radiance  estimates  in  comprehensive  methods  for  thermal- 
structure  resolution.  The  reliabilities  can  also  be  used  to  select  spots 
having  highest  resolution,  for  retrievals  of  individual  thermal  profiles. 

A  major  benefit  of  a  scanning  sensor  is  its  resolution  of  gradient 
along  scans.  The  new  procedure  addresses  this  information  component  of 
the  VTPR  system  in  producing  the  full  scan  for  each  channel. 


3.2  The  Two-State  Interpretation  Model 

An  interpretation  model  forms  the  basis  for  diagnosing  the  VTPR 
radiances  for  their  clear-column  components.  The  simplest  such  model 
is  the  two-state  model. 


The  two-state  model  allows  any  spot  to  be  a  mixture  of  clear 
portions  and  of  cloudy  portions  for  which  all  cloud  top  is  at  the  same  level. 
Furthermore,  the  effectiveness  of  the  model  depends  on  having  pairs  of 
spots,  in  proximity,  for  which  the  cloud  top  is  at  the  same  level  but  which 
differ  in  cloud  proportion.  These  and  other  conditions  are  involved  in  the 
effectiveness . 


Based  on  this  model  an  arbitrary  VTPR  radiance  may  be  expressed  as 
a  combination  of  a  cloud -top  component  and  a  clear-column  component: 
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where  S’ 

subscript  v 
subscript  n 
x 

superscript  P 
superscript  c 
superscript  * 


denotes  the  radiance  expressed  as  a  normalized  radiance 
anomaly, 

denotes  an  arbitrary  channel, 
denotes  an  arbitrary  spot, 
denotes  the  proportion  of  cloudy  area, 
denotes  a  VTPR  radiance, 

denotes  the  cloud-top  radiance  component,  and 
denotes  the  clear-column  radiance  component. 


Expressed  in  terms  of  cloud  proportion,  Eq.  (5)  takes  the  form 


x 


n 


T  *  -  tp 

P,n _ v ,  n 

IT*  -  ZC 
v,  n  v ,n 


(6) 


The  cloud  proportion  applies  to  every  channel  of  the  n-th  spot. 


It  is  convenient  to  introduce  a  vector  notation:  Any  parameter 

which  takes  on  a  value  for  each  channel  of  a  spot  may  be  expressed  as 

P 

a  vector  for  that  spot — one  element  per  channel.  For  example,  ^  n 
denotes  the  set  of  VTPR  radiances  for  the  n-th  spot: 


(7) 


Accordingly  we  define  the  vectors  O  and  0  : 

~  n  ~n 


<p  s  ir  *  -  <i r 

v,  n  v,n  v,  n 


(8) 


e  s  T*  -  Tc 
v ,  r.  v ,  n  v ,  n 


(9) 
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It  is  apparent  from  Eq.  (6)  that  these  vectors  are  parallel: 


<P 


n 


x 


e 

n  ~  n 


(10) 


According  to  Eq.  (5)  the  VTPR  radiance  difference  between  two 
spots,  n  and  n*-l  ,  is  expressed  by 


r  -  3 1  =  x  + 

l\n+l  v  ,n  n+1  p,n+l 


(*  “  Xn+l)  ''  i/,n+l 
"  \K,n  -  ~  \)K,n  * 


(ID 


This  equation  may  be  rewritten  in  the  form: 


v  ,n+l 


-  3' 


V  ,n 


=  (x  -  x  ..  )  (%*  -  \ 

\  n  n+1 J  ^  v,n  v ,n ) 


+  x  I  Tc  -  Tc 

n+1  \  i^,n+l  u,n 


(‘  '  xn+l)  fan*!  -  ^.n) 


(12) 


Of  the  three  terms  on  the  right-hand  side  of  Eq.  12  the  first  term  will 
dominate  if 


(1)  there  is  a  significant  difference  in  cloud  amount  from  one  spot  to  the 
other, 

(2)  there  is  a  significant  difference  between  the  clear  and  the  cloud-top 
radiance  components,  and 

(3)  there  is  relatively  little  change  in  each  of  the  two  radiance  components 
from  one  spot  to  the  other.  This  implies  that  the  cloud  tops  are  at  the 
same  level  in  both  spots. 
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The  probability  that  these  conditions  are  satisfied  is  enhanced  if  the  two 
spots  are  in  proximity.  Except  that  in  regions  of  extensive  overcast  the 
first  condition  is  not  satisfied.  For  those  channels  for  which  a  significant 
portion  of  the  clear-column  component  is  contributed  by  the  underlying 
surface,  the  probability  is  enhanced  over  ocean  areas. 

The  larger  the  VTPR  radiance  difference  between  spots  the  more 
likely  it  is  that  the  right-hand  side  is  dominated  by  the  first  term. 

Based  on  these  conditions,  Eg.  (12)  reduces  to 


y,n+l 


x  -  x  ,  )  6 
n  n+l  /  v,  n 


Introduce  another  vector,  I p  ,  with  clem  -nts 

~  n 


v,  n+l 


Equation  (13)  states  that 


l (i  =  ( x  -  x  ^  0 

~  n  \  n  n+l  J  n 


The  vector  ip  is  parallel  to  the  vector  0  n  ,  in  the  same  or  in  opposite 
direction. 


Equations  (10)  and  (15)  combine  to  yield 


n  i 

-  x  .  ~  n 
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The  vector  p  is  the  adjustment  to  be  added  to  the  VTPR  radiance  vector 
~  n 

for  the  n-th  spot,  to  form  the  clear-column  radiance  vector: 


3*  =  3 P  +  p 

—  n  ~  n  ~  n 


(17) 


The  vector  lb  is  an  ambient  spot-to-spot  difference  between  the  measured 
~  n 

VTPR  radiances.  Under  suitable  conditions  the  vector  ib  gives  an 

■*-  n 

approximation  of  the  vector  axis  for  the  adjustment  vector  p  ^  .  The  scalar 
proportionality  factor  of  Eq .  (16)  is  not  directly  resolved  by  the  VTPR 
radiances . 

In  the  ambience  of  the  n-th  spot  several  estimates  of  the  adjustment 
vector  axis  may  be  formed  using  the  spots  from  a  single  scan: 

From  single  differences: 


P  P 

3  -3 

~  n +m  —  n+m  - 1 


for  m  =  -1 ,  0 ,  1 ,  2 


From  double  differences: 


rr  P  -  3 P 

~  n+m  ~  n+m  -2 


for  m  =  0,1,2 


From  triple  differences: 


3  P  -  3  P 

~  n  im  ~  n+m-3 


for  m  =  1,2 


The  differences  need  not  necessarily  involve  the  n-th  spot.  Spots  of 
adjoining  scans,  if  in  suitable  proximity,  may  be  used  to  form  additional 
estimates  of  the  adjustment  vector  axis. 
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In  order  to  combine  several  estimates  of  the  adjustment  vector 
axis  to  form  a  single  estimate  for  the  n-th  spot  it  is  convenient  to 
normalize  each  estimate  to  unit  length  and  to  set  ito  sign  so  that  it  has 
a  positive  component  along  the  positive  orientation  of  the  vector  space. 
We  define  the  positive  orientation  of  the  vector  space  by  a  vector  which 
has  +1  in  each  element: 

I  =  (1,  1.  ...  1) 


The  normalized,  sign-adjusted  vector  is  denoted  by  N  : 
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~  n 
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(18) 


A  preliminary  resultant  can  be  formed,  from  several  individual 
estimates  of  the  adjustment  vector  axis,  by  weighted  combination.  The 
magnitude  of  the  difference  vector  has  some  value  as  a  weight  because 
it  is  indicative  of  the  desired  interpretation  conditions.  We  have  chosen 
to  weight  by  magnitude  squared.  This  formulation  is  given  by 


N 

~n 


(19) 


The  left-hand  side  is  the  resultant  of  the  several  estimates  combined  on  the 
right-hand  side  by  weighting  according  to  magnitude  squared.  Many 
refinements  can  be  introduced  to  form  the  resultant  by  processes  of  selecting 
and  weighting  individual  estimates.  These  processes  are  aiso  included  in 
the  reevaluation  operation. 


-48- 


The  unit  vector  N  is  an  estimate  of  the  axis  of  correction  from 
~  n  - 

measured  cloud-contaminated  to  clcar-column  radiances  at  spot  n  .  The 

components  are  denoted  by  N^  ^  .  The  amount  of  displacement  along  this 

correction  axis;  i.e.  ,  the  length  of  the  correction  vector,  for  spot  n  ,  is 

not  directly  revealed.  We  denote  this  length,  which  may  turn  out  to  be 

positive  or  negative,  by  L  . 

n 

The  interpretation  model  yields  the  following  information  statements 
at  each  spot,  n  : 


=  L  N  ,  for  all  V  .  (20) 

v ,n  n  v,n  ' 

The  exploitation  of  these  information  statements  is  discussed  in  the  following 
section. 


3 . 3  Formulation  of  the  Error  Functional 

Program  CLRX  is  a  new  application  of  our  Fields  by  Information 
Blending  (FIB)  methodology.*  This  information  blending  methodology,  in 
general  concept,  combines  independent  weighted  estimates  into  the 
non-independent  resultants  implied  by  the  focus  of  all  input  information. 
It  handles  linear  interrelationships  and  other  extensions.  The  present 
application  involves  a  new  extension. 


Holl,  Manfred  M.  and  Bruce  R.  Mendenhall,  1971;  "Fields  by  Information 
Blending,  Sea-Level  Pressure  Version" ,  Final  Report,  Project  M-167, 
Contract  No.  N663 14-70-C-522 6  (Fleet  Numerical  Weather  Central), 
Meteorology  International  Incorporated,  Monterey,  California,  66  pp.  plus 
Appendix. 
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The  blending  of  all  independent  elements  of  information  is 
accomplished  by  minimization  of  a  total  error  functional.  This  error 
functional  is  the  total  weighted  sum  of  the  squares  of  all  resultant 
disparities  with  the  information  estimates.  Each  weight  is  the  inverse  of 
the  error  variance,  or  sum  of  contributing  error  variances,  associated  with 
the  particular  information  estimate. 

The  present  version  of  CLRX  is  based  on  the  following  formulation 
of  the  error4  functional: 


A  FIB  convention  in  notation  is  to  denote  the  resultants  of  the 
blending  by  superscript  asterisks.  The  above  formulation  of  the  error 
functional  is  expressed  in  terms  of  the  resultant  correction  elements, 

p 

O*  ,  to  be  added  to  the  measured  normalized  radiances,  T  ,  to 
V,n  '  v,n 

form  the  resultant  clear-column  normalized  radiance  anomalies: 


T  *  s  UP  f  <p* 

v, n  v,n  P,n 


(22) 
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The  above  error  functional  includes  provision  for  the  following 
elements  of  information: 

(a)  Direct  estimates  of  correction  elements: 

<P  represents  an  assembled  estimate  for  channel  v  and  spot  n  ,  and 
v ,  n 

^  represents  the  assembled  weight  for  that  estimate. 

In  the  present  formulation  of  CLRX  non-zero  weights  occur  only  in 
conjunction  with  zero  estimates  of  the  correction  element.  These 
occurrences  are  assignments  based  on  quantitative  measures  of  the 
smoothness  in  the  radiances  in  the  ambience  of  the  element.  The 
assigned  weight  is  inversely  proportional,  within  bounds,  to  radiance 
variabilities  in  the  ambience.  Low  variability  would ,  generally,  be 
the  case  in  a  clear  atmosphere.  For  channels  with  a  significant  ground 
contribution  the  clear-column-smoothness  expectation  is  generally 
limited  to  oceanic  regions  .  Low  variability  can  also  be  the  case  over 
a  uniform  solid  cloud  top.  Any  assigned  weight  is  tentative  and  all 
weights  are  subject  to  two  cycles  of  reevaluation. 

(b)  First-difference  estimates: 

u  represents  an  assembled  estimate  of  the  first  difference: 

v ,  n 

(p*  -  (0* 

v ,  n+ 1  v ,  n 

^  represents  the  assembled  weight  for  that  estimate. 

In  the  present  formulation  of  CLRX  we  impose  the  expectation  of  low 

variability  on  the  clear-colurnn  radiances: 

*1.  n  “  ^,n+l  '  for  any  u  and  n  .  (23) 
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This  assumed  condition  is  also  implicit  in  how  we  have  formulated  (c) 
below.  The  imposition  of  this  expectation  on  channels  with  a 
significant  ground  contribution  limits  the  applicability  of  the  present 
CLRX  formulation  to  oceanic  regions.  The  imposition  of  Eq.  (23)  on 
the  correction  elements  related  by  Eq.  (22)  takes  the  form  of  the 
estimates , 


v ,  n 


= 


ar 


p 

P,n+1 


(24) 


The  assigned  weight  is  proportioned  to  estimations  of  the  variability 
in  Eq.  (23).  In  the  present  formulation  a  uniform  weight  is  assigned. 


A  refinement  is  introduced  in  the  second  and  third  cycles  of  blending. 
A  reduced  average  value  of  the  clear-column  radiance  gradient  which 
emerges  from  the  first  cycle  is  introduced  into  the  second-cycle 
assembly,  and  that  from  the  second  into  the  third.  This  correction 
takes  the  form  of  a  factor  added  to  the  Eq.  (24)  estimate: 
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(25) 


The  assigned  weight  is  increased  in  proportion  to  the  expected  associated 
reduction  in  variance.  Refinements  such  as  this  one  are  easily 
formulated  and  are  as  easily  improvable. 

(c)  Correction  vector  estimates: 

N 

v ,n  represents  an  assembled  estimate  of  the  component  for 

channel  V  of  the  unit-vector  correction  axis  for  spot  n  . 

C  W  represents  the  weight  appropriate  to  the  product  Ln  N„  „  . 
n  V  ,n  11  l/«n 
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The  basis  and  derivation  of  correction  axes  has  been  discussed  in 
Section  3.2  where  the  information  element  is  expressed  by  Eq.  (20). 
Additional  details  of  the  present  formulation  of  CLRX  are  given  in  the 
Appendix. 

The  resolution  of  the  value,  for  each  spot  n  ,  is  expressed  in 
the  following  section. 


3.4  3.4  The  Blending  System  of  Equations 

The  blending  equations  are  based  on  the  minimization  of  the  total 
error  functional,  E  ,  expressed  by  Eq.  (21).  They  are  obtained  by  setting 


?>E 


b<p* 

v.n 


0 


(26) 


for  each  v  and  n  .  This  procedure  produces  the  general  equation 


A  ( to*  -  ip  \  -  B  ( to*  -to* 

v,n\  v,  n  v,nj  t'.nV  t>,n+l  u,n 

+  -  K, n-1  '  %,n-l) 

+  C  W  ( <P*  -LN  )  =  0 

n  v  ,n  \  v ,n  n  v , n  / 


(27) 


Since  the  Ln  values  are  completely  unspecified  we  are  free  to  use 
them  to  also  minimize  the  error  functional.  By  sotting 
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for  each  n  ,  we  obtain  the  general  equation 


W. 


v  ,ri 


L 
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(29) 


Equation  (29)  defines  each  Ln  in  terms  of  the  resultants  for  that  column: 
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N  W  O* 
y,n  v,n  v  ,n 


where 
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n 
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(30) 


(31) 


The  blending  system  of  equations,  represented  by  Eq.  (27),  may  also 
be  expressed  with  Ln  eliminated,  as  represented  by  Eq.  (30).  The  result 
is 


A ,  (<P*  -  <p  ) 

v ,n  V  u,n  v ,n/ 

+  B  ( O*  -  (p*  ,+W  ) 

v,n  \  v,r\  t»,n+-l  u,n/ 

+  Bi/,n-l  (^t',n  ^ v,n-l  ~  ^  v ,n-l  ) 

+  C  W  ( (0*  -  N  H  X  N  W  o*  \ 

n  t^,n  V  U  ,n  u,  n  n  uu,n  u>,n  tu,n/ 


=  0 


(32) 


in  which  u)  has  been  introduced  as  a  dummy  variable  for  v . 


The  most  economical  solution  of  the  blending  system  appears  to 
be  an  iterative  method  of  successive  approximation:  successive  over¬ 
relaxation  (SOR)  by  columns  (i.e.  ,  spots),  with  the  inner  solution  for 
each  column  by  explicit  inversion.  The  scheme  we  have  devised  uses 
Eqs.  (27)  and  (30)  rather  than  (32).  The  convergence  rate  can  be  judged 
on  successive  values  of  Ln  at  each  n  . 

The  system  of  blending  equations,  as  represented  by  Eq.  (32), 
may  be  expressed  in  matrix  notation: 

MO*  =  F  (33) 

where  O*  is  the  desired  resultant  expressed  as  a  vector  containing  the 
v,n  elements  in  some  orering,  M  is  the  square  coefficient  matrix,  and  F 
is  the  forcing  vector.  In  general,  the  FIB  methodology  produces  coefficient 
matrices  which  are  symmetric  and  positive  definite. 

The  FIB  methodology  also  provides  for  the  resultant  resolution 
weights,  which  measure  the  focus  of  all  information  on  each  solution 
element.  The  resultant  weights  are  represented  by 

A  to  be  associated  with  <P 
V  ,n  t>,n 

According  to  the  FIB  methodology*  the  diagonal  elements  of  the  inverse 
matrix 


contain  the  individual  elements 


★ 

Ibid,  page  49. 
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The  inversion  of  large  matrices  is  not  always  practical  in  a  routine 
operation  stressing  economy.  In  such  applications  TIB  uses  a  perturbation, 

•Up 

or  leverage,  approach  to  establish  estimates  of  the  elements  of  A  .  In 
the  present  situation  the  number  of  elements  are  reasonable  and  the  form 
of  the  matrix  M  is  simple.  Explicit  calculation  of  the  diagonal  elements 
of  the  inverse  matrix  is  warranted. 

Additional  details  including  reevaluation  formulations  are  given  in 
the  Appendix. 


3 . 5  A  Note  about  the  Resultant  Reliabilities 

The  input  to  CLRX  is  in  the  form  of  a  full  scan  of  normalized  radiance 

anomaly  values.  The  weights  assigned  to  the  inherent  information  elements 

are  determined  by  specified  factors  and  diagnostic  schemes  internal  to  the 

program.  These  are  subject  to  adjustments  and  modifications — a  tuning 

process.  The  objective  is  to  have  A*  n  be  a  realistic  weight  for  the 

resultant  correction  element,  <p*  ,  and  equivalently,  by  Eq.  (22),  for 

v ,  n 

the  resultant  clear-column  normalized  radiance  anomaly,  IT 

v  ,n 

We  must  emphasize  that  A*  represents  the  reliability  of  the 

v  ,n 

resultant  normalized  radiance  anomaly.  The  variance  associated  with  that 
normalized  radiance  anomaly  is  given  by 


The  clear-column  radiance  anomaly,  e*  r  ,  is  related  to  the 

normalized  clear -column  radiance  anomaly.  S’*  ,  by  a  normalizing  factor: 

v,n 


e*  =  f  3* 
v,n  v,n  v,  n 


(34) 
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where  ^  denotes  the  normalizing  factor  for  channel  V  for  the  nadir 

angle  of  spot  n  .  The  corresponding  reliability  of  c*  is 

v  ,n 


f 


-2 

v.n 


(35) 


And  if  the  radiance  anomaly  is  to  be  interpreted  as  an  estimate  of 
a  linear  combination  of  thermal-structure  parameters  such  as  expressed 
by  Eq.  (1) — based  on  physical  relationships  and  modelling  approximations  — 
then  the  weight  must  be  reduced  by  the  addition  of  the  contributing  variances 
inherent  in  these  approximations: 


Associated  Weight 


(36) 


where  o  represents  the  sum  of  the  appropriate  estimates  of  contributing 
variances. 


3 . 6  Results  of  CLRX 

The  R&D  version  of  CLRX  includes  options  for  diagnostic  printout 
which  permits  following  of  each  stage  of  the  calculations  through  the 
three  cycles.  This  output  is  too  voluminous  for  inclusion  here. 

The  input  to  CLRX  is  in  the  form  of  a  full  scan  of  normalized 

P 

radiance  anomalies,  H  ,  such  as  depicted  in  Figs.  10-12.  The 

v  ,n 

present  version  of  CLRX  doss  not  include  provisions  for  exploitation  of 

additional  independent  radiance  estimates. 


^577 
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The  standard  output  includes  the  resultant  fields  of  the  dear- 

column  normalized  radiance  anomalies ,  ,  and  the  associated 

v,r\ 

reliabilities.  A*  .  The  input  and  output  fields  can  be  had  in  the  form 
v ,  n 

of  Varian  Plots  (Figs.  14-16)  and/or  in  the  form  of  printouts  (Tables  5-14). 

In  order  to  simplify  the  boundary  treatments,  two  spots  on  each  end  of 
the  scan  are  omitted  from  the  blending. 

Figure  13  shows  the  location  of  the  series  of  scans  which  are 
here  reproduced  as  examples.  The  variability  in  yield  amongst  these 
samples  is  obvious.  The  results  must  speak  for  themselves.  We  have 
not  entered  into  any  evaluation  studies  embracing  correlations,  tunings, 
calibrations,  quality  analyses,  etc.  Evaluation  should  be  performed  under 
real,  or  simulated,  operational  conditions  in  which  all  available  relevant 
information  is  exploited. 

The  present  version  of  CLRX  has  not  been  optimized  for  running 
time.  Such  optimizations  will  be  dependent  on  refinements  and  extensions 
of  the  formulations.  The  approximate  running  time,  in  terms  of  CP  time  of 
a  CDC-6500  computer,  for  the  processing  of  one  scan  through  three  blending 
cycles  and  an  A*  calculation  is  now  33  seconds.  The  expensive  portion 
is  the  A*  calculation  which  uses  about  2  6  seconds;  we  believe  this 
calculation  will  be  amenable  to  considerable  reduction. 
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Fig.  13  Locations  of  the  Consecutive  Complete  Scans  Numbered  18  through  2  7 
(March  26,  1974) 
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Fig.  15  The  Diagnosed  Clear-Column  Radiance  Anomalies ,  T  ,  corresponding  to  Fig.  14.  Two  spots 
at  the  beginning  and  two  at  the  end  are  omitted  in  the  blending. 


19.6  18.2 


Fig.  16  The  Diagnosed  Reliability  Weights,  A  ,  Associated  with  the  Corresponding  Elements  of  Fig.  15 


Tabic  5 


RFC  « ll  NO,  ia 

THb  VTPR  N'ORMALIZI-U  RAUIaMCF  ANUMALILS  (CLRX  INPUT) 


spo  r 

NO, 

CHANNL.  NO. 

1 

2 

3 

4 

5 

6 

W  I  NUUW 

1 

5. 75 

1.55 

-.64 

.73 

7.11 

-  ?  ,  to  to 

*5,21 

2 

a .  2  5 

1.29 

-1.11 

.63 

8,89 

3  ,4o 

7,51) 

3 

6.0  2 

1,0  4 

-2.  U4 

.10 

9.07 

,83 

6,4b 

4 

5,70 

.61 

-1.63 

-.28 

8.4? 

,46 

7,11 

5 

5,33 

.79 

-1.87 

,0b 

8.  U6 

.23 

6,7  8 

6 

6.15 

.93 

-1.4? 

-.05 

8.30 

.  3  6 

8,16 

7 

6,31 

.49 

-2.26 

-.5/ 

8.14 

.56 

8,b6 

3 

5.81 

-.08 

-2.41 

-1.00 

5.76 

-2,44 

2 ,89 

9 

5,23 

.64 

-2.5? 

-.72 

4.47 

-4,21 

-  ,  6  U 

10 

5,35 

.03 

-2.62 

-1,5/ 

.92 

-8,92 

-11,53 

11 

5,40 

.73 

-2,69 

-1.16 

1.68 

-7,to5 

*9,09 

12 

5,43 

.08 

-2.7? 

-.76 

5,38 

-2,64 

2,23 

13 

5,44 

.09 

-2.74 

-.30 

7.22 

-.37 

to,  /8 

14 

5.4? 

,07 

-2.05 

-.71 

7,75 

.4.3 

b.  16 

15 

6,64 

,06 

-2.67 

“1,11 

7.41 

-.25 

7,11 

1  6 

5,3? 

.67 

-2.5B 

-.85 

6,75 

-  .  8l 

3 ,  b  2 

17 

6,49 

.62 

-3.13 

-1.18 

7 . 98 

.49 

10,27 

18 

6.33 

.54 

-2.34 

-1,30 

7.45 

.  0  2 

V.tol 

19 

6,85 

.44 

-2,83 

-.32 

6.50 

-1,20 

3,29 

20 

6,67 

,  86 

-2.63 

-1,34 

6,19 

-1 . 15 

3,95 

?1 

7.0  7 

1 . 36 

-2.97 

-1.66 

0 . 99 

-  .23 

8,56 

22 

6,80 

1.71 

-2.73 

-1.98 

6.53 

-1,09 

6,7  8 

23 

7.69 

2.76 

-3.13 

-2.10 

5 .74 

-1,49 

to,  78 

24 

7,87 

2,95 

-2.87 

-  ?  ,  1 8 

4.68 

-  ?  ,  4  3 

0,7  8 
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Table  5  (Continued) 


REC  >HD  NO,  18 

THE  DIAGNOSED  CLE AR-COLUHN  COMPONENTS  (CL.HX  OUTPUT) 


SPIT 

NO, 

CHANNE.  NO, 

1 

2 

3 

4 

5 

6 

W J  NDUW 

1 

0,00 

o.ou 

o .  on 

n .  ou 

U  ,00 

n ,  oo 

0,00 

2 

0  ,  UO 

o  .no 

0.00 

0,00 

o .  on 

0,00 

0,00 

3 

5,79 

.*2 

-1.82 

-.09 

8.67 

.5/ 

...  7,00 

4 

5,81 

.72 

-1.75 

-.17 

8,46 

,49 

7,12 

5 

5.88 

.76 

-1.78 

-.16 

8.23 

.  39 

7 , 23 

6 

6,04 

.75 

-1.74 

-.27 

8,0  3 

.2/ 

7  ,  o5 

7 

0,17 

,4« 

-2.25 

-  ,56 

7.85 

,1/ 

7,45 

8 

5,97 

.31 

-2.33 

-.69 

7,77 

.12 

7,44 

9 

5,72 

.32 

-2.43 

-.73 

7.69 

,0/ 

7,43 

10 

5,57 

.24 

-2.54 

-.81 

7.63 

,  0  3 

7,4  0 

11 

5,50 

.20 

-2.6? 

-.03 

7.56 

-,01 

7,38 

12 

5,45 

.09 

-2,67 

”,  78 

7,4  9 

-  ,  05 

7,35 

13 

5,45 

.07 

-2,70 

-.79 

7,45 

-.09 

7,31 

14 

5.72 

.11 

-2.43 

-  ,  8  2 

7,4  4 

-  ,11 

7,24 

15 

6,22 

.12 

-2,61 

-1,06 

7,40 

-,2l 

7,21 

16 

6,25 

.39 

-2.65 

-1.04 

7,31 

-.26 

7,22 

17 

6,35 

.53 

-2.72 

-1 , 10 

7.25 

-  ,  3U 

7.21 

18 

6,46 

.57 

-2,57 

-1,19 

7.11 

-.41 

7.21 

19 

6,65 

,68 

-2,69 

-1,25 

6,96 

-.52 

7,20 

20 

6,75 

.93 

-2.73 

-1,44 

6,83 

-.62 

7,19 

21 

6.90 

1.15 

-2.86 

-1,62 

6.78 

-  ,  6t> 

7.19 

22 

6,85 

1.3l 

-2.77 

-1,86 

6,71 

-.79 

7.19 

23 

0,00 

0,00 

o.un 

0 ,00 

0.00 

0,00 

o,  on 

24 

0,00 

0,00 

o .  uo 

0  ,00 

0.00 

0,00 

0,00 

THE 

ASSOCIATED 

RELIABILITIES  (CLHX  OUTPUT) 

1 

0,00 

0.0  0 

u .  uo 

0,00 

0.00 

0 ,00 

U  ,  00 

2 

0.00 

o.oo 

u.  uo 

0 .00 

0  .Of) 

o.uo 

O.UO 

3 

44,89 

60 , 31 

24.95 

41,81 

28.55 

29,10 

15,02 

4 

47.03 

72.0  4 

40.34 

54,82 

36.99 

41.0o 

19,4? 

5 

38.45 

66,66 

48.49 

53,24 

41.49 

45,60 

17,38 

6 

58,87 

110.65 

86.35 

63,43 

42,21 

36,54 

15,4? 

7 

138,46 

115,07 

145,69 

84 , 60 

32.98 

25,68 

12,30 

8 

53,60 

44.22 

55,53 

38.74 

18.91 

16,25 

9.54 

9 

37,73 

30.86 

43,82 

28,47 

14.74 

13,09 

a.ia 

J  0 

35,81 

26.72 

35.4  6 

24,59 

13.59 

12.15 

7 ,34 

11 

41,50 

29,54 

37.82 

29,4/ 

14.61 

12,79 

7,02 

12 

70.80 

55,39 

60.10 

44,78 

19,08 

15,60 

6 , 99 

13 

374,91 

564 ,06 

306,85 

170,86 

35.67 

23,95 

7,22 

14 

59,33 

121,41 

78.40 

104,95 

39,53 

27,69 

6,83 

15 

53,20 

126.33 

81.91 

147,20 

60*12 

39,95 

6,28 

16 

31,63 

50,76 

4  4.02 

51,59 

30 . 21 

25,01 

5.80 

17 

34.50 

7  4,07 

44.48 

52,11 

27.04 

22.53 

5.50 

18 

49,65 

116,18 

55.69 

55,03 

28.49 

23,06 

5,33 

19 

53,33 

47.88 

48.61 

39,38 

25.99 

21 . 4V 

5.20 

20 

54.83 

39.09 

59.15 

41,69 

27,07 

22,46 

5,17 

21 

79,79 

35,  vy 

141.82 

52,73 

31,20 

23,65 

5,22 

22 

51.65 

24,32 

184.20 

56,16 

23.97 

19,91 

5.25 

23 

0,00 

0,00 

0.00 

0,00 

0 .00 

0 ,00 

u.uo 

24 

0.00 

0  ,ou 

O.UO 

0  ,uu 

o.oo 

0 , 00 

0,00 

Table  6 


RECORD  MO.  19 

THE  V  T  P  P  NORMALIZED  RADIANUE  ANOMALIES  (CLHX  INPUT ) 


SPOT  C  H  A  N  N  b  _  MO. 


NO  • 

1 

2 

3 

4 

5 

6 

W  INUUW 

1 

7.29 

1,44 

-1.07 

1 . 56 

9,41 

,/u 

7,4  4 

2 

6,65 

1 .18 

- .  78 

1.35 

8.62 

*  2  4 

6,88 

3 

6.33 

.93 

-1,05 

1.42 

8.80 

1 , 28 

8,49 

4 

6,63 

.50 

-1.30 

1,03 

9,13 

1,6  7 

9.15 

5 

6,26 

.68 

-1.54 

,88 

8.76 

1,51 

0,82 

6 

5.8  3 

.16 

-1.74 

.74 

9 .00 

1,40 

9,64 

7 

6.00 

.27 

-1,93 

.81 

8.83 

1.76 

9,87 

8 

6.1? 

-.19 

-1.42 

.38 

8.87 

l,6l 

10,60 

9 

5,5  9 

-.13 

-1.54 

,55 

8.53 

1,59 

1  U  .  2  0 

10 

o,29 

-.08 

-1.64 

.78 

8,68 

1.7/ 

10,20 

11 

5.7? 

-.06 

-1.71 

.11 

9.00 

2,0  2 

10,93 

12 

6.37 

-.03 

-1.74 

.02 

9.33 

2,33 

11.65 

13 

5.75 

-.02 

-1.76 

-  «  0  2 

9.29 

2,67 

11,65 

14 

6,98 

-.04 

-1.73 

.  55 

9,38 

2.  /3 

11,98 

15 

6,33 

-.06 

-1.69 

.16 

9.05 

2,42 

11,66 

16 

6.89 

.46 

-1.59 

.42 

8.82 

2,61 

12,31 

17 

6,80 

•  4  U 

-2.16 

-  ,  4  U 

9.10 

2,42 

12,31 

18 

7.31 

.32 

-2.66 

-.51 

8,58 

2,33 

12,31 

19 

7,1? 

.80 

-2.39 

-.63 

8.15 

1,94 

11,26 

20 

6,35 

.76 

-2.19 

-1,04 

7.42 

1 , 26 

9,87 

21 

7.3  3 

1.26 

-2.64 

-.96 

7.25 

,22 

9,64 

22 

7.11 

1,60 

-2.39 

-1.1/ 

6,27 

-  ,  63 

6,38 

23 

7.3  3 

1,99 

-2.80 

-2.4/ 

4,04 

-2,94 

4  ,  bl 

24 

8.17 

2.86 

-3.09 

-3.01 

2.49 

-5 ,06 

-  ,  66 
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Table  6  (Continued) 


RECORD  NO,  19 

THE  DIAGNOSED  CLEAR-COLUMN  COMPONENTS  (CLRX  OUTPUT) 


SPOT 

NO, 

channel  no, 
1 

2 

3 

4 

5 

6 

WINDUW 

1 

0,00 

0,  OU 

o,  on 

0.  00 

0,00 

0  ,  UU 

0,00 

2 

0,00 

0,00 

u.un 

0,00 

0.00 

0,00 

0,00 

3 

6,32 

.76 

-1.22 

1.2  7 

8,92 

1,55 

9,03 

4 

6,41 

,S7 

-1,33 

1,07 

9.01 

1,58 

9.15 

5 

6,24 

.48 

-1.52 

.91 

8.91 

1.55 

9,31 

6 

6,04 

.27 

-1.70 

.78 

8,94 

1,53 

9,60 

7 

6,02 

.16 

-1,80 

,73 

8.87 

1,70 

9,91 

8 

6,03 

-.02 

-1.64 

.57 

8.83 

1,73 

10,19 

9 

6,03 

-.07 

-1,61 

.51 

8.82 

1,78 

10,40 

10 

6,12 

• ,  0  6 

-1.67 

.48 

8.90 

1,92 

10,65 

11 

5,99 

-  .  04 

-1.71 

.19 

9,03 

2,07 

10,99 

12 

6,09 

- ,  0  3 

-1.73 

,12 

9.19 

2,27 

11,33 

13 

6.14 

-.02 

-1,76 

.13 

9,26 

2,50 

11,59 

14 

6,34 

-.02 

-1.76 

.20 

9,22 

2,55 

11,74 

15 

6,51 

.07 

-1,73 

,13 

9,05 

2,49 

11,85 

16 

6,69 

.26 

-1.81 

,00 

8.90 

2,49 

12,03 

17 

6,87 

.36 

-2,13 

-.29 

8.82 

2,38 

12,12 

18 

7.12 

.  4  U 

-2.41 

-  ,  4  8 

8.58 

2,31 

12,15 

19 

7.12 

,63 

-2,39 

-.53 

8,46 

2,26 

12,10 

20 

7,00 

.61 

-2,31 

-.70 

8.33 

2,23 

12,17 

21 

7,06 

.77 

-2.41 

*  ,  72 

8,31 

2,22 

12,26 

22 

7,04 

.87 

-2.35 

- »  74 

8,30 

2,25 

12,38 

23 

0,00 

0,00 

0.00 

0,00 

0,00 

0,00 

U,U0 

24 

0,00 

0,00 

0 .00 

0  ,  00 

0.00 

0,00 

u,uo 

THE  ASSOCIATED  REL  1  Afc)  1 1_  I  T  I  ES  (CLRX  OUTPUT) 


1 

0,00 

0.00 

0.0  0 

0,00 

0.00 

0,0U 

0,00 

2 

0,00 

0.00 

0.01) 

0,00 

0.00 

0,00 

u.uo 

3 

66,86 

104,83 

79,78 

78,33 

54,40 

30,16 

13,91 

4 

57,84 

87.29 

69,69 

85,62 

74.09 

45 ,87 

21,37 

5 

43,45 

49,52 

56.24 

60,15 

58,04 

56,0b 

25,23 

6 

52.02 

47,89 

88.51 

66,07 

79,91 

62,7  7 

27,48 

7 

60,76 

51.54 

85,50 

71,85 

72.27 

6  3 , 7 1 

28,57 

8 

39,39 

51.04 

51,70 

46,61 

48.66 

51,46 

27,91 

9 

33,63 

62,51 

56.30 

46,46 

44,77 

52.04 

31.87 

10 

41.83 

287.36 

128,62 

63,53 

60*47 

58,27 

31,88 

11 

43,36 

802 , 85 

478.6? 

71,42 

68.64 

60,06 

32,53 

12 

31,4  6 

759,84 

520 .83 

51 ,15 

63,18 

4  9,92 

33,22 

13 

27,39 

>95,87 

600,56 

58,47 

67,62 

47,87 

40,91 

14 

26,65 

64,02 

218.09 

56,33 

48,45 

48,44 

39,44 

15 

31,59 

57,65 

93,41 

48,55 

51.05 

58,68 

38,49 

16 

46,01 

57,73 

5d .  73 

38.81 

50.91 

62,40 

37,07 

17 

53.36 

56,43 

55.85 

45,44 

43.79 

57,73 

33,26 

18 

61,24 

100.58 

84.61 

87,48 

49,15 

56,43 

25,36 

19 

47.06 

77,43 

112.97 

101,50 

38,15 

32.52 

16.88 

20 

32,94 

55,18 

72.71 

53,20 

24,33 

18,61 

11,05 

21 

31,16 

50,38 

58.55 

33,86 

15.77 

12,12 

8,11 

22 

27,09 

31.53 

38,32 

18.69 

10.31 

8,55 

6,24 

23 

0.00 

0 ,00 

0.00 

0,00 

O.UO 

0,00 

0,00 

24 

0.00 

0,00 

0.00 

0,00 

0,00 

0,00 

0,00 
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Table  7 


RECORD  MO,  20 

THE  VTPH  NORMALIZED  HAulANCE  ANOMALIES  (CLKX  INPUT) 


SPOT 

NO, 

CHANNE.  NO. 

1 

2 

3 

4 

5 

6 

WINDUP 

1 

6,97 

1,33 

-.81 

1,96 

9.96 

1,87 

8,29 

2 

5,63 

.*3 

-1.45 

1,76 

10.25 

1 , 7b 

9 ,  U  2 

3 

6,02 

.82 

-1,05 

1,31 

9,87 

2,05 

9,02 

4 

6,32 

.39 

-1.30 

1.43 

9.75 

2,05 

8,29 

5 

6,57 

-  .09 

-1.54 

1,28 

9,62 

1.82 

10,07 

6 

5.52 

.05 

-1.74 

1.14 

9,5? 

2,53 

10,80 

7 

6,31 

.16 

-1.38 

.61 

9,44 

2 , 5  u 

9, 74 

8 

6,44 

-.30 

-1.53 

.77 

9.91 

2 , 5o 

11,13 

9 

5,90 

-,?4 

-1.65 

.95 

10.00 

2 , 7 1> 

11,13 

10 

5,99 

-.19 

-1.75 

,68 

9.71 

2,69 

10,80 

11 

6.0  3 

-  ,  16 

-1.82 

.50 

9.51 

2, 76 

10,80 

12 

5,43 

-.80 

-1.85 

.41 

9.94 

3.07 

11,85 

13 

6.07 

-.79 

-1.67 

.37 

9.89 

3,04 

11 , 45 

14 

5.42 

-.81 

-1.84 

.46 

9.98 

2. 73 

11.13 

15 

6,64 

-  ,  16 

-2.35 

.55 

9.57 

2,79 

11,45 

16 

6,57 

-.87 

-2.25 

.32 

9.42 

2.96 

11.45 

17 

6.49 

-.26 

-2.15 

-.0  0 

9.71 

2,42 

11,85 

18 

5.75 

-.34 

-2.01 

-.11 

9,18 

2,70 

11 , 45 

19 

6,86 

.77 

-1.84 

-.23 

8.24 

1.94 

10,80 

20 

6,67 

.  64 

-2,96 

-1,24 

7.94 

1,18 

10.07 

21 

7,07 

1.14 

-2.75 

-1.06 

6,81 

-.23 

8,o9 

22 

8,04 

1.49 

-2,51 

-1.38 

6.34 

-.72 

7,96 

23 

7.03 

2.43 

-2.25 

-1,5b 

6.00 

-,6b 

8,29 

24 

8,47 

2,74 

-2.54 

-2.61 

4.48 

-2,76 

5,86 

Table  7  (Continued) 


RECORD  NO,  20 

THE  DIAGNOSED  CLEAk-COLUMN  COMPONENTS  (CLHX  OUTPUT) 


SPOT 

NO, 

CHANNEL  NO. 

1 

2 

3 

4 

5 

6 

WINDUW 

1 

0,00 

0,00 

0.0  0 

0,0U 

0.00 

0,00 

0  ,  uo 

2 

0,00 

0  .  ou 

0,00 

0,00 

0,00 

o,ou 

0,0  0 

3 

6,04 

,68 

-1,13 

1,38 

9,89 

2,04 

9,28 

4 

6,15 

-1.30 

1,37 

9.86 

2,01 

9,48 

5 

6,46 

,10 

-1.52 

1,28 

9.80 

2,01 

9,89 

6 

6,14 

.10 

-1,67 

1.11 

9.67 

2,31 

10,23 

7 

6,20 

-.04 

-1.51 

,85 

9,68 

2,4/ 

10 , 43 

0 

6,20 

-,1b 

-1.52 

,78 

9,76 

2,56 

10,65 

9 

6,05 

'.lb 

-1,63 

.80 

9,81 

2,68 

10,83 

10 

5,99 

-.23 

-1,73 

,68 

9,75 

2,79 

10,90 

11 

5,95 

-.32 

-1.81 

.55 

9,69 

2.83 

11,04 

12 

5,91 

» ,  43 

-1.82 

,49 

9,75 

2,89 

11,20 

13 

6,00 

-,53 

-1.8* 

.42 

9,79 

2,90 

11,28 

14 

6.10 

-.51 

-2.04 

,47 

9,76 

2,82 

11,30 

15 

6,26 

-.49 

-2,17 

,44 

9,64 

2,80 

11,37 

16 

6,43 

-,5b 

-2.23 

,26 

9,50 

2.7/ 

11,45 

17 

6,44 

-.40 

-2.17 

,02 

9,46 

2,60 

11,48 

18 

6,12 

-.13 

-2.1? 

-.22 

9.05 

2,48 

11,32 

19 

6,49 

.43 

-1.99 

-  ,  30 

8.76 

2,3l 

11,20 

20 

6,40 

.31 

-2.55 

-.72 

8,58 

2,0/ 

11,09 

21 

6,49 

,36 

-2,90 

-,90 

8,47 

1,95 

11,05 

22 

6,54 

,36 

-2.8ft 

-1,01 

8,42 

1.91 

11, 08 

23 

0.00 

0,00 

0 ,00 

0,00 

0,00 

o.oo 

u,uo 

24 

0,00 

0.00 

0 ,00 

0,00 

0  ,U0 

0,00 

0,00 

THE  ASSOCIATED  RELIABILITIES  (CLHX  OUTPUT) 


1 

0,00 

0 ,  ou 

0,00 

0,00 

0.00 

u.ou 

U  ,00 

2 

0.00 

0  .00 

0.00 

0.00 

0 .00 

0,00 

0,0  0 

3 

38,45 

44,95 

49.83 

46,32 

45.68 

51.50 

IB.  00 

4 

55,79 

54.78 

62.12 

55.11 

57,40 

58,80 

31.2# 

5 

167,07 

82,38 

186,11 

65,59 

84.73 

63,56 

23,83 

6 

48,07 

92,21 

156.92 

50,51 

55.75 

50.87 

21,68 

7 

39,84 

82,50 

96,36 

52.62 

55.35 

53,92 

20,88 

8 

42,44 

104,44 

88.9ft 

74.89 

62.36 

58,26 

24,01 

9 

44,94 

104,87 

69.02 

70,14 

56.35 

53,12 

30 , 44 

10 

45,16 

64,51 

8  7.99 

54.76 

51.58 

52,09 

35,  ia 

11 

48,31 

51,45 

170,44 

61,87 

59,60 

60.29 

32,09 

12 

3  4,7? 

4  3.27 

337,33 

9  4,07 

57.08 

57,73 

30,29 

13 

28.96 

38,45 

74,94 

128,75 

51.39 

52,15 

33,82 

14 

28,31 

44,21 

75,68 

108,06 

59,22 

53,63 

41,00 

15 

31,17 

43,51 

7  7.15 

66,4/ 

59,01 

54,33 

46,40 

16 

44,81 

40.70 

77,57 

59,39 

58,70 

49,64 

46 ,  U# 

17 

173,25 

46.94 

134.76 

55,21 

79,26 

48.39 

42,19 

18 

88,85 

86.12 

158,07 

79.85 

55,04 

43,51 

38,62 

19 

79,54 

102,78 

77,34 

55,17 

40*33 

32,05 

27,43 

20 

58,39 

79,06 

51.52 

44,82 

29.19 

22,08 

17,89 

21 

24,26 

27,44 

5  7,85 

36,03 

16.67 

13,74 

11.64 

22 

14,05 

15.85 

30,52 

21.38 

10,98 

9,48 

8,28 

23 

0,00 

0,00 

0,00 

0.00 

0,00 

0,00 

0,00 

24 

0,00 

0,00 

0,00 

0,00 

0.00 

0,00 

0,00 
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Table  8 

wee  )R I)  MJ,  21 

THfc  VTPP  NU?MALIZf:J  H  A  L)  I  A  N  C  F  ANUrlALItS  (CLKX  INPUT; 


SPOT 

NO, 

CHANNfc. 

1 

NO. 

2 

3 

4 

5 

0 

U  I  N  1 )  0  W 

1 

7.19 

1.22 

-2.62 

.11 

8,40 

15 

4.21 

2 

6.46 

.42 

-2.33 

.01 

8.71 

.63 

5,26 

3 

6,33 

,  06 

-2.60 

-.51 

7,37 

-1.15 

2,83 

4 

6,01 

-.37 

-3.40 

-2.48 

2,40 

-6,50 

-3,50 

5 

5.63 

-.20 

-4.29 

-3,70 

1.56 

-7,06 

-'2,11 

6 

5.83 

-  .t'X 

-3,83 

-2.72 

3.76 

-4,07 

1.77 

7 

6.62 

-.60 

-2.80 

-  .  0  7 

7.97 

.12 

6,  /I 

3 

6.50 

-1.07 

-2.30 

.  09 

8 ,45 

.  0  2 

7.7  7 

9 

5.5? 

-1.01 

-3.0  7 

-1,31 

6.20 

-1,91 

4.21 

10 

5,66 

-.96 

-3,17 

-1,0/ 

7,30 

-  .  oO 

4,61 

11 

6,97 

-.93 

-3,24 

-.67 

8.05 

.46 

7,04 

12 

6.9? 

-.91 

-2.61 

.81 

9.85 

3.0  7 

10,86 

13 

6.33 

-1.45 

-2,63 

.1/ 

10  .  75 

3.41 

11,92 

14 

6.9  3 

-.92 

-2.60 

.20 

1U  .  41 

3,10 

11.92 

15 

6,96 

-.93 

-2.56 

.35 

9,99 

3  ,  03 

11.92 

3.6 

6.83 

-.32 

-2.4  7 

.61 

11.23 

4,03 

14,02 

17 

6.80 

-.37 

-3.03 

.38 

11.52 

4.65 

13.69 

13 

7,32 

-.45 

-2.89 

.28 

11.00 

4 , 19 

13,69 

19 

7,69 

.  ou 

-2 .72 

.  IP 

11,11 

4,05 

14,  U2 

20 

7,51 

.64 

-2.53 

-,25 

10.30 

3,43 

13, 36 

23 

7.23 

1.03 

-2.97 

-,67 

9.  72 

2,86 

12,31 

22 

8.25 

.34 

-2.73 

-.99 

8.7  3 

1 . 93 

11.2« 

23 

7,90 

1.78 

-3.03 

-1,68 

6 . 90 

-1.51 

4,61 

24 

8,0  3 

1.99 

-3,43 

-2.21 

4 ,46 

-3.23 

2,50 
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Table  8  (Continued) 


RFC  >RD  NO,  21 

THE  DIAGNOSED  CLEAN-COLUMN  COMPONENTS  (CLHX  OUTPUT) 


SPOT 

channel 

NO, 

NO, 

1 

2 

3 

4 

5 

6 

WINDOW 

1 

0,00 

0,00 

0,00 

0,00 

0.00 

0,00 

0,  UO 

2 

0.00 

0  ,oo 

0,00 

0,00 

0,00 

o ,  uu 

0,00 

3 

6.41 

.10 

-2.10 

1,00 

9,79 

1,96 

8,75 

4 

6,32 

'.12 

-2,12 

,99 

9,82 

2,06 

9,03 

5 

6,26 

",  32 

-2.17 

.9/ 

9,84 

2,18 

9,33 

6 

6.26 

",  54 

-2.20 

.95 

9,88 

2,28 

9,60 

7 

6,32 

-,69 

-2.20 

,96 

9,91 

2,39 

9,87 

8 

6.22 

",90 

-2,10 

,91 

9,94 

2,50 

10,14 

9 

6,32 

.,93 

-2.36 

,  86 

9,90 

2.67 

10,44 

10 

6,51 

-.94 

-2.58 

.79 

10,05 

2,86 

10,74 

11 

6,73 

.,  96 

-2,68 

.74 

10,10 

3,03 

11,04 

12 

6,85 

-.96 

-2.59 

,74 

10.15 

3,20 

11,32 

13 

6,55 

-1 , 3l 

-2.67 

,  36 

10.44 

3,28 

11,66 

14 

6,87 

-.98 

-2.59 

,34 

10.46 

3,40 

12,05 

15 

6,93 

«.  80 

-2.58 

,40 

10.44 

3,74 

12,42 

16 

6,94 

-.65 

-2.58 

,45 

10.64 

3,9/ 

12,82 

17 

7,01 

-.57 

-2, an 

,34 

10.86 

4,12 

13,19 

18 

7,30 

* ,  30 

-2.81 

,22 

10,90 

4,14 

13,55 

19 

7,53 

.10 

-2,71 

,06 

10,90 

4,2* 

13,80 

20 

7,50 

.32 

-2,61 

-,04 

10,82 

4,15 

14,01 

21 

7,50 

.5-4 

-2.77 

-.2/ 

10,81 

4,23 

14,28 

22 

7,73 

.51 

-2.69 

-,41 

10,84 

4,3* 

14,58 

23 

0,00 

0.00 

0,00 

0,00 

0.00 

0  ,00 

0  ,U0 

24 

0,00 

0,00 

o.on 

0  ,00 

0,00 

0,00 

0,00 

THE 

ASSOCIATED 

RELIABILITIES  (CLHX  OUTPUT) 

1 

0,00 

0,00 

0.00 

0,00 

0.00 

0.0  0 

O.UO 

2 

0,00 

0.00 

0.00 

0,00 

0,00 

0,00 

0,00 

3 

44,15-  . 

45,03 

20.99 

7,64 

4,32 

3,8  7 

3,11 

4 

32,67 

34.72 

17,46 

8,30 

4.83 

4,33 

3,44 

5 

26,92 

33.77 

17.54 

9,41 

5.52 

4,93 

3.85 

6 

26,72 

37,43 

22.21 

11,55 

6.52 

5,78 

4,40 

7 

30,19 

51.29 

43.41 

lo,  08 

8.0  3 

7,04 

5,14 

8 

27.47 

50.26 

51.32 

20,23 

9.69 

8,41 

5,98 

9 

23.19 

41.60 

30,63 

17,60 

10.56 

9,26 

6,65 

10 

23,21 

48.40 

32.56 

19,13 

12.61 

11. OU 

'7,76 

4  4 

L  1 

31.0  7 

54,62 

49,35 

25,91 

17.39 

14  ,  75 

9,68 

12 

92,6? 

255.32 

337.82 

67,26 

34,25 

25,93 

13,61 

13 

140,47 

166,69 

345.14 

66,02 

44,76 

36,71 

19,71 

14 

97,16 

70.98 

671.9* 

69,44 

51.68 

37,83 

26,21 

15 

179,76 

120,13 

116.54 

137,53 

41.91 

34,76 

33,53 

16 

88,33 

115,23 

66.82 

83,24 

34,56 

30,93 

22,86 

17 

63,80 

143,23 

71.06 

94,64 

39,19 

37,68 

29,31 

18 

63,73 

63,62 

99.29 

110,44 

51.81 

43,42 

41,28 

19 

75,50 

61.57 

92.19 

101,14 

47,88 

35,92 

31,45 

20 

61,27 

51,19 

100.65 

75,14 

32.39 

24,72 

20,16 

21 

41,17 

51,94 

83,69 

42,03 

18,88 

14,78 

11.88 

22 

25,90 

38,59 

52.30 

24,62 

11.83 

9,91 

8,24 

23 

0,00 

0.00 

U.00 

0,00 

0.00 

0,00 

0,00 

24 

0.00 

o.ou 

0.00 

0,00 

0,00 

0,00 

0.00 

-70- 


Table  9 


RPCJHH  NO,  22 


THfc  VTPR  NORMAL  IZt:L  KAUIANUE  ANUMALJtS  (CLKX  INPUT) 


SPOT 

C  H  A  N  N  fc  _ 

NO. 

NO, 

1 

2 

3 

1 

7.19 

.79 

3T 

X 

H 

1 

2 

0,56 

-2.11 

3 

6,94 

.2/ 

-1.71 

4 

6.63 

.50 

-1.96 

5 

7.50 

.02 

-1.64 

6 

6,45 

-.39 

-2.40 

7 

7.24 

-,?8 

-2.03 

8 

6.7  5 

-.85 

-2.1'> 

9 

5,59 

-1.34 

-2.30 

10 

6,2? 

-.74 

-2.95 

11 

5.72 

-1.20 

-3.0? 

12 

5,7  4 

-1.24 

-2.51 

13 

6. 33 

-.68 

-1.87 

14 

6,93 

-.70 

-1.84 

15 

6,33 

-.71 

-2.46 

16 

6.83 

- .  7o 

-2.36 

17 

6,80 

-.15 

-2.81 

18 

6,69 

- .  23 

-2.07 

19 

7.17 

.22 

-2,50 

20 

7.61 

.76 

-2.97 

21 

0.77 

.60 

-3.96 

22 

6,4? 

.*5 

-3.73 

23 

6,7  7 

2.  CO 

-2.92 

24 

8,70 

2.21 

-3.21 

4 

5 

0 

k  I  NUUK 

1.89 

12,09 

4,14 

111,  73 

1.76 

11.33 

2, do 

9,28 

1.83 

11.50 

3,05 

8,23 

1,44 

11,35 

3,11 

10.  73 

1,29 

10.45 

2,50 

8,02 

1.1* 

11.10 

3,13 

10,34 

1.11 

10.49 

2,73 

8,95 

1.27 

10.96 

3 , 10 

11.39 

1.44 

10.61 

2,50 

9,08 

-.39 

8,94 

1,20 

8.02 

2.72 

2.11 

-6,44 

-1.52 

.42 

9.59 

1  ,dl 

1U  .  01 

.80 

11.43 

3 . 6  o 

12.44 

.94 

11.01 

4. 06 

12,44 

.55 

9.74 

2,27 

10,73 

.81 

10.97 

3,57 

12,44 

.48 

10.74 

4,13 

12,77 

.37 

10.74 

3 , 3  0 

12.44 

,  30 

10.76 

3,28 

12,44 

-.15 

9.51 

2,23 

11,06 

4,88 

-2.74 

-12, 7  0 

-11,33 

5,73 

-8.26 

-19,87 

=24,96 

2.69 

.68 

-8,62 

*8,17 

2,12 

3,08 

-5,64 

-2,57 

-71 


Table  9  (Continued) 


RECORD  NO,  22 

THE  DIAGNOSED  CLEAH-COLDMN  CUMPUnENTS  (CLRX  OUTPUT) 


SPOT  CHANNEL  NO, 


NO, 

1 

2 

3 

4 

5 

6 

WINDOW 

1 

0,00 

0,00 

0,00 

0,00 

0.00 

0,00 

o ,  un 

2 

0,00 

0,00 

0,00 

0,00 

0.00 

0,00 

0 , 00 

3 

6,64 

.37 

-1.88 

1,61 

11.41 

3,09 

10,22 

4 

6,65 

.31 

-1.94 

1 , 46 

11.29 

3,07 

10,34 

5 

6,57 

,03 

-2.07 

1,32 

11.18 

3,06 

10,45 

6 

6,45 

-  ,  3u 

-2,27 

1.19 

11.10 

3 , 1U 

10,53 

7 

6,36 

-,5  6 

-2,24 

1.1« 

10.99 

3 ,  U  7 

10,68 

8 

6,19 

- ,  8y 

-2.25 

1,28 

10.95 

3,06 

10,81 

9 

5,79 

-1,25 

-2.29 

1,45 

11,02 

3,09 

10,88 

10 

6,02 

-.98 

-2,55 

1,16 

10,99 

3,15 

11,11 

11 

5,97 

-1,05 

-2,46 

1,03 

10.95 

3,18 

11,32 

12 

5,99 

-1.07 

-2.30 

,89 

10.9  0 

3,20 

11,52 

13 

6,32 

-.74 

-1,97 

,73 

10.87 

3,22 

11,68 

14 

6,61 

-.75 

-2,05 

,79 

10.72 

3,32 

11,82 

15 

6,67 

-,69 

-2,25 

,73 

10.71 

3,35 

11,98 

16 

6,79 

-.69 

-2.41 

.72 

10.72 

3 , 3b 

12,13 

17 

6,69 

- ,  32 

-2,68 

.50 

10,62 

3,44 

12,27 

18 

6,77 

- ,  15 

-2,66 

,36 

10,66 

3,29 

12, 3« 

19 

7,17 

.23 

-2,52 

,33 

10,68 

3,26 

12,44 

20 

7,59 

.67 

-2.87 

.21 

10.58 

3,29 

12.53 

21 

7,49 

.'0 

-3.01 

,14 

10.54 

3,31 

12,69 

22 

7,47 

.76 

-3.09 

,08 

10.51 

3,33 

12,80 

23 

0,00 

0,00 

u.  on 

0,00 

0,00 

0,00 

0,00 

24 

0,00 

0,00 

0,0  0 

0,00 

0,00 

0,00 

0,00 

THE 

ASSOCIATED 

RELIABILITIES  (CLRX  OUTPUT) 

1 

0,00 

0,00 

0,00 

0 ,00 

0.00 

0,00 

0,00 

2 

0,00 

o.oo 

0,00 

0,00 

0.00 

0  ,ou 

0 ,00 

3 

24,87 

46,10 

44,39 

33,26 

26,15 

30 . 70 

5,19 

4 

36,92 

58.35 

74.59 

70,39 

35,76 

45.99 

6,41 

5 

28,99 

37,08 

42.52 

54,50 

30.46 

34,  86 

7,60 

6 

39,89 

43,24 

64,80 

123,98 

43.41 

49,47 

9,37 

7 

30,56 

38,64 

53.05 

69,45 

33.18 

34,80 

10,04 

8 

40,77 

58,47 

90.70 

52,47 

34.62 

33,70 

11.6U 

9 

115,33 

231,91 

248,20 

58,16 

26,75 

24,92 

12,46 

10 

52,34 

64.95 

65,77 

30 ,03 

18.55 

17,26 

11,23 

11 

41,06 

45.19 

40 ,68 

25,84 

1 6  .  o  3 

15,28 

11,03 

12 

83,33 

114,39 

90.8^ 

42,06 

19,52 

16,88 

12,1? 

13 

158,06 

295,59 

256,58 

92,93 

26,07 

20,78 

14,48 

14 

62,87 

98.01 

105,06 

109,20 

28,15 

21,75 

16,45 

15 

59,27 

74,85 

53.20 

68,24 

27.50 

21.74 

17,74 

16 

177.29 

110,86 

85.  >  0 

122,62 

42,51 

28,74 

23,23 

17 

178,68 

92,76 

125.90 

114,11 

52.40 

32,69 

28,27 

18 

166,94 

128,70 

220,22 

155,02 

61.72 

39,98 

33,61 

19 

553,50 

535,03 

795.18 

158,41 

44.45 

32,07 

20,60 

20 

175,99 

100.28 

162,64 

44,59 

20.34 

16,42 

11.95 

21 

27,00 

32,28 

26.03 

16,19 

11.25 

9,93 

8,09 

22 

14,79 

16,76 

14.77 

9,87 

7.77 

7,11 

6,11 

23 

0,00 

0,00 

0.00 

0,00 

0,00 

0,00 

0,00 

24 

0,00 

0.00 

0.00 

0,00 

0.00 

0,00 

0,00 

Table  10 

RfcC  i MU  NU.  23 

T  Kit  VTP?  NORMAL  I  ZF.O  RAUIANUE  ANOMALIES  (CLKX  INPUT) 


SPOT  CHANNE »  NO, 


NO, 

1 

2 

3 

4 

5 

6 

WINDUH 

1 

6,37 

.68 

-1.95 

1 . 68 

9,b8 

1,73 

8,23 

2 

6.87 

.53 

-2.2? 

2,06 

11.60 

3,24 

10,34 

3 

6,63 

.17 

-2.48 

2.14 

10.69 

2,6/ 

8,6? 

4 

4,46 

-,«1 

-5.27 

-8.52 

-14.17 

-25,13 

=> 3U  ,23 

5 

5,94 

-.09 

-2.31 

2.09 

9.65 

1,30 

7,90 

6 

6.14 

-  .  50 

-1.96 

It  94 

8.9? 

.81 

6,5? 

7 

5,6* 

-1.05 

-2.14 

.91 

7.36 

-1.08 

3 . 68 

S 

5.81 

-.96 

-2.30 

.48 

5.93 

-2.58 

2 , 0  0 

9 

4.0  3 

-1.45 

-9.19 

-17,19 

-26.03 

-36,93 

-42.48 

10 

4.10 

-1 . 4u 

-9.29 

-17.41 

-27.70 

-38.54 

-45,25 

11 

4.13 

-1.37 

-4.98 

-8.19 

-12.26 

-22,38 

•24,31 

12 

6 ,  U  3 

-.80 

-2,61 

,61 

8.82 

.  62 

b,  b? 

1.3 

6,07 

-1.34 

-2.63 

1,15 

9.29 

1.4l 

V.61 

14 

6.04 

-.81 

-2.60 

1,24 

10*33 

2.95 

11,7? 

15 

6.02 

-  •  82 

-2.56 

.74 

9.48 

2.27 

10.66 

16 

6,57 

-.87 

-2.47 

.51 

9.76 

2.48 

10,66 

17 

o,49 

-.26 

-2.37 

,78 

10.05 

2,64 

10,66 

18 

5.76 

-1.00 

-4.119 

-3,57 

-5,67 

-16, 1U 

•19,76 

19 

4,99 

-  ,  43 

-3.93 

-3.11 

-6,15 

-17,00 

-21,80 

20 

5.43 

-.01 

-3.74 

-4.14 

-8,47 

-19,64 

*26,41 

21 

7,05 

.49 

-2.86 

-1,38 

4,76 

-3,79 

1.25 

22 

6.81 

.  y5 

-3.29 

- ,  60 

7.21 

-.51 

7.1 7 

23 

8.2? 

1.24 

-3.03 

-  .  6  / 

7.86 

-  ,  0  6 

7,90 

24 

7.7  8 

2.75 

-3,32 

-1,81 

4  » 9(1 

-3.32 

3.68 

73- 


Table  10  (Continued) 


NEC  )RD  NO,  23 

THE  0 1 AGNOSED  CLt AH-COlUMN  CUMRUNENTS  (CLKX  OUTPUT) 


SPOT 

CHANNEu 

NO. 

NO, 

1 

2 

3 

4 

5 

6 

WINDUW 

1 

0,00 

o.  no 

o,  un 

0  ,  OU 

0,00 

o,  nu 

0,00 

2 

0,00 

0,  ou 

U  .  00 

0,  OU 

0.00 

0,00 

0 ,00 

3 

6,62 

.16 

-2,48 

2,17 

10.67 

2,60 

8,69 

4 

6,33 

.03 

-2.35 

2,2/ 

10,56 

2,48 

8,84 

5 

6,04 

-.10 

-2,18 

2,38 

10.45 

2,36 

8 , 9H 

6 

6,15 

-.47 

-1.85 

2,37 

10.33 

2,29 

9,  U8 

7 

5.93 

-.79 

-1.75 

2,27 

10,25 

2,24 

9,22 

8 

6,00 

-.89 

-1. 76 

2.12 

10.16 

2,18 

9,38 

9 

6,03 

*-.98 

-1.9? 

1.93 

10.09 

2,15 

9,57 

10 

6,05 

-1,02 

-2.U3 

1.73 

10,02 

2,11 

9,75 

11 

6,09 

-  ,98 

-2.25 

1.53 

9,94 

2,0  7 

9,93 

12 

6,12 

-.89 

-2.43 

1,33 

9,87 

2,03 

1U.11 

13 

6,11 

-1.30 

-2.58 

1,28 

y.75 

2,01 

10,26 

14 

6,12 

•1,04 

-2,60 

1,06 

9,68 

2,10 

10,38 

15 

6,24 

■» ,  9y 

-2.56 

,78 

9,57 

2,22 

1U.48 

16 

6,55 

-,8b 

-2.49 

.54 

9,63 

2,32 

10,49 

17 

6,49 

-.28 

-2.40 

.66 

9,74 

2,32 

10,43 

18 

6,59 

-.15 

-2.59 

,48 

9,66 

2,29 

10,52 

19 

6,75 

.05 

-2.78 

,29 

9.57 

2,25 

10,61 

20 

6,91 

.31 

-2.91 

,10 

9,49 

2,21 

10,70 

21 

7,06 

.61 

-2.98 

-  ,  08 

9.41 

2,18 

10,79 

22 

7,03 

.86 

-3.15 

-.24 

9,34 

2,15 

iu,yo 

23 

0,00 

0,00 

o.uo 

0 ,00 

0.00 

0,00 

0,00 

24 

0,00 

0,00 

u ,  uo 

0,00 

0,00 

0,U0 

0,00 

THE 

ASSOCIATED 

Rfcll  AbILI TIES  (CLRX  OUTPUT) 

1 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

U  ,00 

2 

0.00 

0,00 

0.00 

0,00 

0.00 

0 ,00 

u,uo 

3 

534,55 

815.53 

400.07 

41.39 

8.56 

6,91 

3,69 

4 

46,10 

46,96 

44,16 

23,58 

8.11 

6,75 

3,85 

5 

232.01 

235,36 

172.06 

31,82 

8,87 

7,26 

4,16 

6 

139,39 

117.6Q 

111.65 

27,88 

8.82 

7,31 

4,39 

7 

63.00 

61.36 

51,62 

1 J .  98 

8.43 

7,18 

4,63 

8 

38,47 

43.20 

32.59 

15,44 

8,28 

7,22 

4,89 

9 

22,23 

27.81 

20.69 

13.67 

8.44 

7,48 

5.2S 

10 

20.80 

28,63 

19.  /4 

14,31 

9,26 

8.22 

5,89 

11 

23.34 

34.27 

26,59 

10,22 

11*19 

9,76 

6 ,87 

12 

145,35 

159.72 

106.32 

35,34 

15.97 

13,15 

8,58 

13 

475,74 

472.97 

445.56 

81,55 

25.11 

18,73 

21,05 

14 

68,36 

61.45 

180 . 48 

47,89 

29,64 

20,38 

23,0« 

15 

62,84 

73.01 

179.97 

54,97 

43,18 

29,01 

27.3* 

16 

809,30 

838,22 

832.01 

188,50 

48.86 

35,9/ 

29,57 

17 

636,82 

668.21 

632.13 

137,02 

31,41 

23,23 

12,68 

18 

32,06 

34,74 

32.33 

25,48 

14.83 

12.57 

8,53 

19 

20,64 

25.88 

23.57 

16,76 

10.14 

8,85 

6,4'/ 

20 

19,30 

25.07 

23.85 

14,51 

7,97 

6,97 

5,25 

21 

24,27 

31.34 

32.8? 

14,92 

6.75 

5,85 

4,44 

22 

29,66 

37,50 

48.40 

16,47 

5,87 

5,05 

3,85 

23 

0 ,00 

0,00 

0.00 

0,00 

0,00 

0.00 

0,00 

24 

0,00 

0,00 

0.00 

0.00 

0,00 

0,00 

0,00 

-74- 


Table  11 


RECORD  NO,  24 

THE  VTPR  NORMALIZED  RADIANCE  ANOMALIES  (CLRX  INPUT) 


SPOT 

NO, 

CHANNE.  NO, 

1 

2 

3 

4 

5 

6 

WINDOW 

1 

9,06 

2,79 

-2.01 

-2,03 

1,92 

-7,45 

-6,92 

2 

8,25 

3.12 

-3,61 

-4.21 

-1.86 

-11,91 

-11,46 

3 

8,59 

2,06 

-2.53 

-1,26 

5.35 

-3,29 

3,22 

4 

8,23 

1.71 

-3,43 

-3,63 

-.85 

-10,4b 

•  7,58 

5 

7.21 

1,20 

-4.30 

-7,56 

-12,18 

-23,64 

-27,86 

6 

6,75 

.78 

-3.17 

-3,83 

-5,57 

-16,27 

-ai,54 

7 

6,29 

.22 

-3,33 

-4,7/ 

-9,34 

-20,00 

-25,76 

8 

5,77 

.29 

-2.92 

-3,55 

-3,98 

-14,18 

-16,01 

9 

5,84 

,34 

-3,02 

-3,80 

-3.80 

-14,34 

-14,95 

10 

5,29 

-.27 

-3,66 

-6,69 

-10.16 

-20,95 

-25,03 

11 

5,94 

- » 26 

-3.15 

-3,65 

-2.67 

-12,25 

-12,12 

12 

5,34 

-.24 

-3,72 

-6,26 

-10,32 

-21,05 

-26,08 

13 

5,32 

-  .  p0 

-5,01 

-8 , 36 

-14,57 

-25,91 

-  3 1 , 6  8 

14 

5,30 

-.27 

-5,53 

-10.34 

-16,82 

-28,16 

-34,12 

15 

5,87 

-.30 

-3,60 

-5,4/ 

-11.49 

-22, /2 

-30 ,30 

16 

5.80 

-.35 

-2.96 

-3 , 66 

-5.98 

-16,85 

—23,32 

17 

5,71 

-.41 

-3,39 

-4,93 

-9,06 

-20,09 

.27,53 

18 

5,59 

-.50 

-3.90 

-6,10 

-11,95 

-23,72 

-31,02 

19 

5,42 

.05 

-4,93 

-8,78 

-15.31 

-26,87 

=33,07 

20 

5.22 

-.09 

-4,18 

-6,61 

-11.38 

-22,21 

=28,92 

21 

6,22 

.40 

-3,29 

-4.31 

-8,16 

-19,84 

=26,81 

22 

6,54 

.72 

-3.U5 

-3,4/ 

-6.25 

-18,48 

=25,03 

23 

6,79 

.43 

-3,45 

-2.99 

-4.09 

-16,09 

=20.55 

24 

6,29 

1.24 

-3,19 

-4,09 

-7.63 

-19,80 

=28,19 

-75- 


Table  11  (Continued) 


RECORD  NO,  24 

THE  DIAGNOSED  cle ar-column  COMPONENTS  (CLRX  UUTPUT) 


SPOT 

NO, 

C H A N N E _  NO, 

1 

2 

3 

4 

5 

6 

K1NDUW 

1 

0,00 

0,0  0 

u.uo 

0.00 

0,00 

0 ,00 

0,00 

2 

0,00 

0  .  ou 

u ,  on 

0 ,00 

0.00 

0  ,  OU 

0,00 

3 

7,60 

1.39 

-3,51 

-4,2/ 

-4.20 

-15,05 

-14,89 

4 

7,47 

1,30 

-3.47 

-4,16 

-4,16 

-14,88 

ol5tll 

5 

7.14 

1,10 

-3.34 

-4,01 

-4,10 

-14,70 

-15,35 

6 

6,75 

.82 

-3.16 

-3,84 

-4,05 

-14,55 

-15,57 

7 

6.27 

.*’1 

-3,06 

-3,73 

-4.03 

-14 , 43 

•15,73 

8 

5,82 

,30 

-2.93 

-3,64 

-4,02 

-14,32 

-15,88 

9 

5,83 

.31 

-3,04 

-3,84 

-4,03 

-14 , 35 

•15,89 

IQ 

5,79 

.01 

-3.15 

-3,84 

-3,98 

-14,08 

• 16 , 16 

11 

5,80 

-.22 

-3.20 

-3,83 

-3,93 

-13,80 

•16,44 

12 

5,69 

- ,  26 

-3,06 

-3,60 

-3,83 

-13,54 

•16,82 

13 

5,65 

-.32 

-2.92 

-3,3/ 

-3,73 

-13,2/ 

•i/,21 

14 

5,70 

-  .  3l 

-2.77 

-3,13 

-3,62 

-13,05 

•17. 59 

15 

5,81 

•  ,  32 

-2,64 

-2,90 

-3,52 

-12,82 

.17,97 

16 

5,84 

-.34 

-2.58 

-2,69 

-3,41 

-12,64 

•18,32 

17 

5,84 

-.32 

-2.59 

-2,54 

-3,33 

-12,40 

•18,64 

18 

5,87 

-.22 

-2.57 

-2,39 

-3,25 

-12,24 

•18,93 

19 

5,93 

-.01 

-2,60 

-2,26 

-3,18 

-12,04 

•19,21 

20 

6,19 

.17 

-2,63 

-2,12 

-3.10 

-11,94 

•19,47 

21 

6,4? 

.39 

-2.70 

-2,00 

-3,02 

-11,81 

=19,74 

22 

6,53 

.52 

-2.74 

-1,89 

-2,95 

-11,72 

•19,97 

23 

o,co 

0 ,00 

0,00 

0,00 

0.00 

0,00 

U  ,00 

24 

0,00 

0.00 

0,00 

0,00 

0.00 

0,00 

0,00 

THE 

ASSOCIATED 

RELIABILITIES  (CLKX 

OUTPUT) 

1 

0,00 

0.00 

0.00 

0,00 

0,00 

0 ,00 

0,00 

2 

0,00 

0,00 

0.00 

0,00 

0.00 

0 ,00 

0 ,00 

3 

14,96 

14,30 

13.17 

5,67 

1.94 

1.71 

.74 

4 

24,1  6 

24,38 

23.09 

7,16 

2.09 

1.83 

.76 

5 

26.54 

29,25 

24,39 

8,67 

2.25 

1.96 

,78 

6 

37,36 

40.92 

39,06 

11,37 

2.43 

2,11 

,80 

7 

42,45 

43,80 

43.64 

15,31 

2.59 

2,27 

.82 

6 

224,33 

433,41 

331.00 

30,49 

2.82 

2.47 

,86 

9 

324,22 

317,07 

25/. 90 

19,04 

2.74 

2,36 

,89 

10 

48,20 

48,32 

42.19 

12,44 

2.56 

2.21 

,9Q 

11 

37,41 

81,99 

50.74 

11,16 

2.45 

2.11 

,91 

12 

38,79 

39,33 

23,94 

8,26 

2.30 

2.00 

,91 

13 

31,70 

29,55 

17.06 

6,74 

2,17 

1,92 

.91 

14 

29,74 

32.19 

15.58 

5,90 

2.08 

1.85 

,92 

15 

36,49 

40,09 

1/.25 

5,43 

2.01 

1,80 

.92 

16 

5*1,97 

61,56 

22,93 

5,14 

1.95 

1.76 

.92 

17 

41,31 

41,49 

19.69 

4,71 

1.86 

1.70 

.91 

18 

31,79 

32,13 

15,78 

4,32 

1,76 

1.65 

.90 

19 

24,63 

29,51 

14,51 

4,02 

1.70 

1.60 

,88 

20 

21,75 

28,33 

15,48 

3,80 

1,64 

1,55 

.87 

21 

23,33 

30,16 

17,67 

3,60 

1.58 

1.51 

,8b 

22 

23,67 

25,06 

17.09 

3,34 

1.51 

1.45 

,88 

23 

0,00 

0.00 

u.uo 

0,00 

0,00 

0,00 

0,00 

24 

0,00 

0,00 

u.un 

0,00 

0.00 

0.00 

0,00 

-76- 


Table  12 


RECuHD  NO,  25 

THE  VTPR  NORMALIZED  RADIANCE  ANOMALIES  (CLRX  INPUT) 


SPOT 

NO, 

CHANNEL.  NO, 

1 

2 

3 

4 

5 

6 

W  I  N DU W 

1 

9,95 

3.33 

-2,12 

-.59 

7,01 

-1.20 

1.31 

2 

9,16 

2,46 

-2,39 

.0/ 

7.57 

-.29 

2,7  0 

3 

8.89 

1.51 

-2. 09 

-  .  35 

7,84 

-.26 

4 , 80 

4 

8,54 

1.71 

-2,33 

-.63 

6,47 

-2,  IU 

,98 

5 

8,14 

.66 

-2.54 

-.20 

5 , 69 

-2,99 

1,64 

6 

7,70 

.78 

-3,28 

-3,43 

-  ,  69 

-10,29 

*9,95 

7 

7,84 

,  88 

-2.90 

-2,31 

2,12 

-6,79 

•  3,63 

8 

7,43 

•  2  V 

-4,23 

-5.32 

-5.45 

-15,45 

=15,94 

9 

6,26 

,34 

-4,33 

-7,63 

-11.82 

-22,89 

=27,20 

io 

6,33 

-.27 

-5.08 

-8.94 

-14,90 

-26,07 

•  32 , 4  7 

u 

5.73 

*  ,  26 

-5.11 

-9,01 

-14,02 

-24,56 

-28,59 

12 

6,39 

- .  7y 

-6,34 

-11 .14 

-18,83 

-30,31 

o  3  6 , 3  6 

13 

5.7  4 

-.25 

-4,47 

-7,97 

-12.68 

-23,46 

-28,26 

14 

6,34 

-.27 

-4,44 

-8.38 

-15,44 

-26.46 

-32,80 

15 

6.29 

-  .30 

-4.3/ 

-7.23 

-10,54 

-21,01 

-25 , 48 

16 

5,59 

-.35 

-4.27 

-6,99 

-9.35 

-20,12 

-23,6b 

17 

6,12 

-.«i 

-4.16 

-7,19 

-13,30 

-24,49 

»3l , U9 

18 

6.62 

-  .5u 

-4. 01 

-5.7U 

-10,56 

-21,18 

-26,87 

19 

5,83 

-.61 

-3,83 

-6. 30 

-11.98 

-23,19 

-29,31 

20 

6,26 

-.09 

-3,63 

-4.61 

-7.95 

-19,35 

-24, U4 

21 

6,64 

.  *0 

-2.74 

-1,80 

-1,78 

-12,65 

=17,00 

22 

5.72 

.17 

-3,82 

-6.42 

-12.15 

-24,50 

=32 , 1 4 

23 

5,96 

.43 

-7.89 

-16,41 

-27.50 

-40,38 

-49,00 

24 

6,60 

.70 

-4.53 

-6,99 

-12.67 

-26,01 

-35,30 

-77- 


Table  12  (Continued) 


RECORD  NO,  ?5 

THE  DIAGNOSED  CLE AH-COLUMN  COMPONENTS  (CLKX  UUTPUT) 
SPOT  CHANNEw  NO. 


NO, 

1 

2 

3 

4 

5 

6 

1 

0,00 

0,  ou 

U,  00 

0,00 

0,00 

0,  uo 

2 

0,00 

0.0U 

0  ,  UO 

0  ,  OU 

0.00 

0,00 

3 

8,16 

1 , 07 

•  3  ,  U  2 

-3,66 

-1.81 

-11,19 

4 

8,03 

.97 

-3.10 

-3,69 

-1.53 

-10,85 

5 

7,85 

.83 

-3,22 

-3,72 

-1.25 

-10,49 

6 

7,65 

.75 

-3.3  8 

-3,7/ 

-1.00 

-10,19 

7 

7.55 

,  66 

-3,48 

-3,85 

-.46 

-9,49 

8 

7,41 

.49 

-3,67 

-3,91 

.06 

-8,83 

9 

7,07 

.35 

-3,77 

-3,97 

.65 

-8,10 

10 

6,76 

.18 

-3,83 

-4,00 

1,16 

-7,48 

11 

6,52 

,  06 

-3,85 

•  4,05 

1,69 

•  6,83 

12 

6,33 

«  ,05 

-3.93 

-4,02 

2.12 

-6,32 

13 

6,27 

*.12 

-4.02 

-4,01 

2.58 

-5,77 

14 

6,25 

-.19 

-4.11 

-3,89 

2.90 

-5,38 

15 

6,22 

-.24 

-4,13 

-3,81 

3.27 

-4,94 

16 

6,13 

-.27 

-4,07 

-  3 , 6  U 

3.49 

-  4  ,  ’  >  7 

17 

6,23 

•  ,  26 

-3.93 

-3,40 

3,76 

-4,35 

18 

6,30 

-.21 

-3,71 

-3,08 

3.90 

-4,17 

19 

6,33 

-.10 

-3,40 

-2,79 

4,08 

-3,95 

20 

6,40 

.05 

-3,04 

-2,41 

4,15 

•  3,86 

21 

6,49 

.20 

-2,65 

-2,09 

4,27 

-  3 , 7  0 

22 

6,32 

.15 

-2,59 

-2,06 

4,56 

-3,36 

23 

0,00 

0  ,  OU 

0,00 

0,00 

0.00 

0 ,00 

24 

0,00 

0.00 

0,00 

0,00 

0.00 

0 ,00 

THE 

ASSOCIATED 

RELI  AHILI TIES  (CLKX 

OUTPUT) 

1 

0,00 

0 ,00 

U  .00 

0 ,00 

0.00 

0 ,00 

2 

0,00 

0,00 

0 .00 

0,00 

0,00 

0,00 

3 

13,91 

12.25 

12,46 

3,46 

1.02 

.81 

4 

21.35 

17.55 

18.07 

3,95 

1.05 

.83 

5 

33,77 

29,7l 

25.28 

4,55 

1.10 

,  86 

6 

121,66 

144,85 

52,68 

5,34 

1.14 

,89 

7 

42,95 

49.80 

30.27 

4,84 

1.11 

,87 

a 

32,20 

39,56 

26,50 

4,44 

1.00 

,86 

9 

23.12 

28.78 

2U  .36 

3.96 

1,06 

.84 

10 

21.53 

24,69 

15,99 

3,60 

1,03 

.83 

u 

19,74 

23.93 

14,17 

3,35 

1, 0 

,82 

12 

2  0,14 

22,58 

14.21 

3,15 

.99 

.81 

13 

20,80 

26,49 

16.84 

3.01 

.93 

.80 

14 

23.59 

30,41 

23,62 

2,90 

.99 

,79 

15 

23,43 

32,09 

29.04 

2.83 

.95 

,78 

16 

21.33 

31,20 

30.95 

2,75 

.94 

.77 

17 

22,30 

29.80 

29.69 

2,69 

.93 

.77 

18 

2  *  ,  4  2 

27,31 

27.36 

2,65 

,92 

.76 

19 

23,06 

26,47 

24.16 

2,63 

.91 

.75 

20 

30.48 

31.27 

21,08 

2.61 

.90 

.75 

21 

42,19 

45,98 

19.21 

2,58 

.90 

,74 

22 

18,76 

23,73 

11.24 

2,35 

.87 

,72 

23 

0,00 

0,00 

0,00 

0,00 

0.00 

0.00 

24 

0,00 

0,00 

0,00 

0,00 

0.00 

0,00 

-78- 


WINDUW 
U,  UO 
0,  UO 
=>11,68 
-12,  2a 
*12,19 
-23, 3D 
•14,29 
•15,25 
•16, 2« 
-r/.iv 

•ia,i3 
-ia,vi 
-IV, 7* 
-30, 3« 
•31, U7 

•  21 , b  0 

•  32,19 
•32,63 
-33,11 
•33, «7 
.33,69 
-24,44 

0  ,  UO 
U  ,  UO 


U  ,  UO 
U  ,  UO 
,43 
.43 
.44 
.  4t> 
.49 
.  4t> 
,45 
.44 
.  4  > 
,44 
,44 
,44 
.44 
.44 
.44 
.44 
.44 
.44 
.44 
.43 
0.00 
U,  UO 


REC.'WD  ?6 


Tabic  1  3 


T.-tb  VTPR  NORMALIZE  KAuIAMCE  ANUMALItS  ( U  L  K  X  INPUT) 


SPOT 

NO, 

CHAN Mb  _  Nu, 

1 

2 

3 

4 

5 

6 

R  I  NDUW 

1 

1  0 . 6  0 

3.86 

-1.35 

,02 

8,17 

,  6  6 

2.17 

9,13 

3.01 

-1.6? 

.04 

9.25 

1.06 

6,38 

3 

0,90 

2.61 

-1.87 

,73 

8,53 

1,02 

3,95 

4 

b ,  5  5 

2,26 

-1.45 

1.06 

9,53 

2.34 

V,  15 

5 

6,77 

1 ,  iu 

-1.67 

-.21 

7 , 78 

.26 

5,66 

6 

7,7  0 

.67 

-3.72 

-7.11 

-6,26 

-15,06 

-16,93 

7 

7.2? 

.77 

-4.54 

-8.62 

-10,56 

-19,92 

-20,82 

8 

6,81 

,1b 

-4.0? 

-  7 , 8d 

-9,95 

-19,7/ 

-2U , 62 

9 

7.41 

.23 

-3.48 

-3.42 

i ,  2  n 

-6,99 

-4,75 

10 

7,43 

,28 

-1.70 

.05 

7.57 

-  ,  36 

3,22 

3.1 

8.13 

-.26 

-2.38 

-.62 

5,59 

-2,09 

-.27 

12 

5,7  6 

-.24 

-4.9  3 

-8.02 

-10.58 

-20 ,56 

-23,25 

13 

5,74 

-.25 

-4.90 

-10,02 

-15,26 

-25,76 

-29,84 

14 

7,49 

-.2/ 

-3,02 

-3.69 

- .  03 

-8,73 

-6 , 85 

15 

7.43 

.25 

-2.29 

-.87 

5.85 

-2,15 

1 , 18 

3  6 

7,36 

.  20 

-2.85 

-1,10 

4,74 

-4,20 

-2,61 

17 

7.2  6 

.14 

-1.42 

1.28 

9,84 

2,05 

6,  U 5 

18 

7,76 

.05 

-1.26 

1.72 

10.33 

2,74 

6 , 71 

19 

7,59 

.60 

-2.40 

-2,01 

3.09 

-4.86 

-3,03 

20 

7,3  9 

,46 

-1.53 

1 .30 

10.19 

1  .  75 

4,61 

21 

8 ,39 

.83 

-1.3o 

1.54 

11.09 

3.8/ 

1U.53 

22 

8,70 

1 . 26 

-1.71 

1,82 

11.75 

4 , 2o 

12,64 

2  3 

6,3  0 

.96 

-1.45 

1.86 

12.15 

4 , 86 

14 , 75 

24 

8,42 

1.77 

-1.18 

1.96 

11.86 

4,96 

14,75 

-79- 


Table  13  (Continued) 

RECOHD  MO,  ?6 

THE  DIAGNOSED  CLb  AK-COlUMN  COMPONENTS  ( CL.KX  OUTPUT) 


SPOT 

NO, 

CHANNEL  MO, 

1 

2 

3 

4 

5 

6 

WINDUW 

1 

0,00 

0,00 

0,00 

0  ,00 

0,00 

0,00 

‘0,00 

2 

0,00 

0,00 

u.un 

0,00 

0,00 

0  ,  ou 

0,00 

3 

8,71 

2,16 

-1.73 

,65 

8,64 

1,28 

6,37 

4 

8,57 

1,94 

-1,61 

,50 

8,64 

1,27 

6,46 

5 

8,76 

1.^4 

-1.54 

.37 

8.72 

1,32 

6,46 

6 

8,45 

.97 

-1,56 

,45 

8,84 

1,39 

6,40 

7 

8,13 

.74 

-1.58 

,54 

8,95 

1,46 

6,33 

8 

7,96 

.51 

-1.59 

,63 

9,06 

1,54 

6,28 

9 

7,81 

.34 

-1,59 

,72 

9,10 

1,61 

6,22 

10 

7,75 

.25 

-1,54 

,82 

9,29 

1,69 

6,19 

11 

7,94 

.07 

-1.64 

,91 

9,39 

1,77 

6,17 

12 

7,85 

.04 

-1.65 

.97 

9,50 

1,85 

6,18 

13 

7,73 

.05 

-1.66 

1,04 

9,60 

1,93 

6,19 

14 

7,69 

.09 

-1,67 

1,10 

9,70 

2,01 

6,21 

15 

7,56 

.18 

-1.67 

1,18 

9,81 

2,10 

6,23 

16 

7.41 

.17 

-1.54 

1.29 

9,91 

2,1b 

6,25 

17 

7,30 

.13 

-1.38 

1,38 

10,01 

2,26 

6,27 

18 

7.72 

,08 

-1.33 

1,53 

10,11 

2,36 

6,27 

19 

7,65 

.32 

-1,39 

1,54 

10,34 

2,43 

6,28 

20 

7,53 

.49 

-1,44 

1,56 

10,57 

2,49 

6,27 

21 

7,91 

.85 

-1.45 

1,51 

10.67 

2,54 

6,35 

22 

7,93 

.71 

-1,55 

1,61 

10,81 

2.6U 

6,36 

23 

0,00 

0,00 

0,0  0 

0,00 

0,00 

0  ,00 

0,00 

24 

0,00 

0 ,00 

0,00 

0,00 

0,00 

0,00 

0 ,00 

THE 

associated 

RELIABILITIES  (CLHX  OUTPUT) 

1 

0,00 

0.00 

0,00 

0 ,00 

0.00 

0 ,00 

0,00 

2 

0.00 

0 .00 

0.00 

o.ou 

0.09 

0,00 

0,00 

3 

65,25 

29,99 

126,56 

24,22 

11.42 

8,25 

2,79 

4 

134,73 

82,27 

132,20 

27,49 

10,47 

7,91 

3,11 

5 

201,05 

156.64 

150.91 

26,76 

9.45 

7,39 

3,29 

6 

32,12 

37.19 

27,75 

15,49 

7.6.) 

6,31 

3,21 

7 

23,51 

30.87 

19.40 

12,33 

6. 7  2 

5,69 

3,1* 

8 

23,25 

30.13 

18.71 

11.49 

6,25 

5,35 

3,10 

9 

31,99 

40.01 

23,99 

12,12 

6,10 

5,20 

3,19 

10 

63,46 

82,81 

50.82 

14,44 

6.18 

5,21 

3,2« 

11 

33,45 

41,77 

31.05 

12,31 

5,90 

5,02 

3,25 

12 

20,89 

26,20 

19.62 

10,66 

5.72 

4,90 

0,2? 

13 

19,23 

24,73 

17.68 

10,32 

5,77 

4,93 

3,33 

14 

24,02 

28,63 

20,47 

11,04 

6,05 

5,10 

3,4s* 

15 

41.14 

44,73 

30.33 

13,09 

6,62 

5,44 

3,60 

16 

51,35 

52,55 

33.83 

16,68 

7. 56 

5 , 96 

3,33 

17 

499,91 

647,91 

140.92 

30,89 

9.42 

6,91 

4,15 

18 

402,41 

448,20 

170.39 

30,18 

9,65 

7,03 

4,05 

19 

54,66 

52,50 

44.55 

21,73 

9,18 

6,46 

3,60 

20 

106,71 

144,33 

154.12 

29,88 

10,14 

6,45 

3,30 

21 

46,35 

67,67 

70 .00 

31,54 

10.31 

5,63 

2,94 

22 

22,12 

26,28 

33,42 

23,60 

8,71 

4,79 

2,62 

23 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

24 

0,00 

0,00 

o ,  no 

0,00 

0.00 

0,00 

0,00 

-80- 


Li  R  D  NO.  2  7 


Table  14 


T+£  VrPR  NORMALIZED  RADIANCE  ANOMALIES  (CLRX  INPUT) 
S»Of  CHANNEL  NO. 


NO. 

1 

2 

3 

4 

5 

6 

WINOON 

9.65 

4.21 

-1.35 

-.13 

3.07 

2.02 

9.41 

i. 

5.55 

3.  e9 

-1.U6 

.45 

9.61 

2.59 

9.08 

> 

9.21 

2.  S4 

-1.57 

.  13 

8.68 

1. 40 

6.65 

u 

9.46 

1.93 

-1.56 

.36 

9.89 

2.27 

9.08 

5 

5.46 

2.00 

-1.7  5 

.29 

9.62 

1.  73 

8.75 

6 

5.01 

1.  55 

-1.97 

-.  87 

6.5  3 

-1.39 

4.14 

7 

5.16 

1.65 

-2.68 

-3.99 

.29 

-7.84 

-4.61 

6.5'- 

1.17 

-4.67 

-8.47 

-9.52 

-19.17 

-19.63 

9 

6.57 

.  56 

-5  .43 

-10.37 

-12.61 

-2  2.37 

-22.79 

1  " 

6.0  2 

.61 

-5.52 

-11.39 

-14.73 

-25.18 

-27.73 

11 

6.67 

.62 

-4.89 

-8.04 

-6.6  9 

-18.33 

-20.03 

12 

5 .70 

•  64 

-3.73 

-7.  05 

-6.88 

-16.08 

-16.47 

X  J 

6.0  5 

.63 

-6  .76 

-13.83 

-13.13 

-29.39 

-32.63 

14 

6.66 

.  61 

-3.02 

-3.  21 

.31 

-8.44 

-6 . 66 

lb 

6.60 

.58 

-2  .9  5 

-3.  02 
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4. 


Retrieval  Demonstration 


« 

n> 


The  present  section  is  included  in  order  to  illustrate  that  we  have 
developed  the  major  technical  components  for  an  operational  capability  to 
exploit  satellite-observed  sounding  data  for  resolution  of  atmospheric 
thermal -structure  variabilities.  We  include  a  demonstration  of  thermal- 
structure  retrieval  using  clear-column  radiances,  diagnosed  by  CLRX, 
as  input  to  the  general  structure-blending  retrieval  capability  developed 
under  Phase  One.* 

This  demonstration  required  modification  of  the  structure-blending 
retrieval  scheme  from  the  13-level,  11-layer  atmospheric-structure  model 
used  in  Phase  One,  to  the  18-level,  16-layer  model  used  in  Phase  Two. 

Two  examples  are  selected:  one  spot  from  Record  No.  18  and  one 
spot  from  Record  No.  27 — both  nearest  to  zero  nadir  angle  for  the  scan. 

The  CLRX-produced  radiances  and  weights  must  be  transformed,  as  shown 
in  Tables  15  and  16,  for  input  to  the  retrieval  scheme.  The  input  first- 
guess  estimates  and  weights  for  the  atmospheric  thermal-structure 
parameters  and  sea-surface  temperatures  are  based  on  information 
independently  provided  to  us . 

The  standard  output  for  these  two  retrieval  examples  are  given  in 
Tables  17  and  18,  respectively.  In  both  examples  channels  722.0  and  746.5 
are  rejected  as  gross  errors  in  the  reevaluations.  These  values  do  not  check 
out  on  the  basis  of  all  information  (i.e. ,  weighted  estimates)  and  the 
relationships.  The  implication  is  that  channels  are  not  reading  as 
advertized  by  their  specified  transmission  functions.  This  is  a  plausible 
conclusion.  Expectations  based  on  Figs.  2  through  9  are  not  bome  out 
by ,  e .  g . ,  Figs .  1 0  and  1 1 . 

★ 

Detailed  and  documented  in  the  Task  Three  Report  referenced  on  page  5. 
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Table  15 


Spot  13  <«l  K'-cord  No.  19.  Columns  (2)  and  (5)  are  :  ' 
CI.KX  products.  Columns  (4)  and  (8)  express  the  inp„ 
radiance  information  to  the  structure-blending  retrieval  sr 


(1) 

(2) 

(3) 

(4) 

(5) 
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49 

2.04 

722.0 

9.26 
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10.77 

67.62 

,49.99 

'  *»4 

1.52 

746.5 

2.50 

1.350 

3.38 

47.87 

26.27 

''  »l 

1.20 

834.0 

11.59 

1.518 

.  17.59 

40.91 

17.75 

0() 

0.96 

■'0  0.25 


I.eqond  to  the  columns: 

(l)  Channel  identification 

(7)  CLRX-diagnosed  clear-column  normalized  radiance  anoma'..  . 

from  Table  6  xtracted 

(9)  The  normalizing  factor  extracted  from  Table  4 

(4)  The  radiance  anomaly  adjusted  according  to  Eq.  (34) 

(-,)  CLRX-diagnosed  reliability  associated  with  column  (2),  e\i 

Table  6  ''‘"ed  from 

((,)  The  reliability  adjusted  according  to  Eq.  (35) 

(/)  Example  of  estimated  variances,  cr2  ,  contributing  to  Eq  .  v -  \ 

explained  in  Section  3.5  '  « 

(i!)  The  reliability,  adjusted  according  to  Eq.  (36),  for  associ.i 

column  (4),  as  input  to  the  retrieval  scheme  '  with 
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Table  16 


Spot  13  of  Record  No.  27.  Columns  (2)  and  (5)  are  the 
CLRX  products .  Columns  (4)  and  (8)  express  the  input 
radiance  information  to  the  structure-blending  retrieval  scheme. 
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Legend  to  the  columns: 

(1)  Channel  identification 

(2)  CLRX-diagnosed  clear-column  normalized  radiance  anomalies  extracted 
from  Table 

(3)  The  normalizing  factor  extracted  from  Table  4 

(4)  The  radiance  anomaly  adjusted  according  to  Eq.  (34) 

(5)  CLRX-diagnosed  reliability  associated  with  column  (2),  extracted  from 
Table  14 

(6)  The  reliability  adjusted  according  to  Eq.  (35) 

o 

(7)  Example  of  estimated  variances ,  a  ,  contributing  to  Eq.  (1)  as 
explained  in  Section  3.5 

(8)  The  reliability,  adjusted  according  to  Eq.  (36),  for  association  with 
column  (4),  as  input  to  the  retrieval  scheme 


-85- 


cycle  number  two  |  Cycle  number  three 
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Cycle  number  two  I  cycle  number  thrfe 


5. 


Design  of  the  Basic  Comprehensive  Capability 


Supplemental  developments  and  refinements  are  necessary  in  order 
to  establish  the  basic  capability  for  processing  satellite  VTPR  sounding 
data  for  atmospheric  thermal-structure  resolution.  The  required  components 
and  interfaces  are  shown  schematically  in  Fig.  17.  Commentary  follows. 
Components  are  identified  by  number  in  the  figure. 

Component  1: 

The  major  requirement  is  for  independent  thermal-structure  estimates, 
to  be  used  in  diagnosing  the  radiances  and  for  purging  the  diagnosed 
correction  vectors  in  CLRX,  and  for  blending  with  the  radiances  in  the 
Structure  Blending  operation  (7).  A  complete  set  of  thermal-structure 
parameters  is  required  for  every  spot  of  a  scan.  The  time  of  the  scan 
determines  the  timeliness  of  the  charts  required.  The  position  of 
successive  spots  are  transformed  into  I,J  locations. 

New  subroutines  developed  by  Mil  for  FNWC  will  be  used  for 
interpolation  to  I,J  locations.  These  are  identified  as  INTRP  and  INTRPS. 
Both  are  based  on  zero-order  and  first-order  spatial  continuity,  INTRPS 
is  for  application  to  ocean  parameters — in  the  present  context,  sea-surface 
temperature — and  involves  a  spatial-covariance-dissociation  (SCD)  field 
demarcating  oceanic  continuity.  The  interpolation  will  also  be  designed 
to  establish  that  the  scan  is  in  an  oceanic  region. 

Transforms  will  be  required  to  transform  the  available  set  of  FNWC 
thermal-structure  parameters,  into  the  set  selected  for  the  exploitation  of 
the  VTPR  system. 

Component  2: 

The  first-guess  thermal-structure  scan  sets  must  be  transformed 
into  first-guess  values  of  normalized  clcar-column  radiance  anomalies. 

The  necessary  transforms  for  a  VTPP  system  of  specified  transmission 
functions,  can  be  produced  by  the  capability  described  under  Section  2.1 
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above,  for  a  selected  thermal-structure  parameterization.  These  transforms 
will  have  to  be  produced  as  required. 


The  addition  of  refinements,  such  as  correcting  the  sea-surface 
component  for  attenuation,  will  depend  on  the  availability  of  such 
techniques . 

Component  3: 

The  scan  radiances  are  transformed  into  normalized  radiance 
anomalies,  using  the  capability  described  under  Section  2.4  above. 

These  transforms  will  have  to  be  produced  as  required. 

Component  4: 

The  basic  capability  for  diagnosing  the  clear-column  radiance 
components  from  cloud-contaminated  VTPR  measured  radiances ,  will  be 
refined  to  exploit  first-guess  radiances  and  the  implied  gradients  across 
the  scan.  The  first-guess  radiances  will  also  be  used  in  the ‘diagnosis, 
weighting,  and  purging,  of  information  elements  in  the  scan,  and  will 
also  serve  as  first-guess  values  to  accelerate  the  convergence  of  blending 
by  iterative  techniques . 

Component  5: 

A  simple  scheme  will  be  designed  for  lateral  checking  of  the 
diagnosed  clear-column  radiances,  exploiting  continuity  with  preceding 
and/ or  following  scans .  This  reevaluation  technique  can  make  use  of 
purported  reliability  weights. 

Component  6: 

A  high-quality  sample  of  diagnosed  clear-column  radiances  will 
be  collected  for  calibration  studies.  The  selections  will  be  based  on 
high  reliability  weights.  Values,  weights,  position  and  time  da.a  will 
be  saved. 
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Component  7: 

Structure  Blending  applies  to  individual  spots  (i.e.  ,  columns)  of 
the  scan.  The  retrieval  technique  is  one  dimensional.  Blending  retrieval 
can  be  limited  to  only  those  spots  having  highest  reliability  in  any  portion 
of  a  scan,  and  relative  to  adjoining  scans,  or  every  spot  of  a  scan  may  be 
retrieved.  The  combination  coefficients,  which  relate  a  radiance  to  the 
thermal-parameter  set,  differs  for  each  spot  dependent  on  the  nadir  angle. 
This  dependency  will  have  to  be  incorporated . 

The  clear-column  radiances  are  diagnosed  and  weighted  by 
component  (4).  The  first-guess  thermal-structure  estimates  produced  by 
component  (1)  must  have  associated  reliability  weights.  Class  weights 
may  be  used.  Refinements  may  be  added  for  estimating  associated  and 
contributing  variances . 

The  combined  resultant  thermal-structure  parameters  and  sea-surface 
temperature  estimates  and  associated  reliability  weights  are  produced.  The 
implication  of  the  independent  radiances  are  also  expressed. 


The  development  of  the  basic  capability  will  involve  evaluations  at 
various  stages  of  component  operations. 


Instead  of  listing  all  the  improvements  and  extensions  that  have 
come  to  mind  for  these  developments — in  this  already  overly  lengthy 
report — we  will  make  do  with  one  sweeping  statement:  The  improvement 
of  applications  of  the  FIB  methodology,  is,  by  the  very  nature  of  the 
methodology,  a  continuous  constructive  process. 

Acknowledgment.  —  A  sense  of  appreciation  behooves  me  to  acknowledge 
Mrs.  B.  A.  Hawkins  for  all  programming,  and  Di .  W.  F.  Weigle  for  writing 
the  Appendix. 
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APPENDIX  A 


The  Component  Operations  of  FIB-CLRX 

1 .  Introduction 


The  general  Fields  by  Information  Blending  (FIB)  methodology  has 
been  developed  by  Mil  and  successfully  used  for  the  analysis  of  horizontal 
distributions  of  both  scalar  (e.g.,  sea-level  pressure,  sea-surface 
temperature)  and  vector  (e.g. ,  surface  wind)  fields.  For  the  current 

application  it  has  been  designed  to  analyze  the  two-dimensional  distribution 

*  P 

of  corrections ,  <P  ,  to  the  measured  radiance  anomalies,  3"  ,  to 

v ,  n  v ,  n 

produce  the  clear-column  anomalies,  ^  ^  n  • 

The  following  sections  will  discuss  the  form  of  the  error  functional 
used  in  this  application,  the  resultant  blending  system  of  equations,  the 
methods  of  solution  of  the  blending  system,  the  assembly  and,  finally,  the 
reevaluation  and  recycling  of  the  analysis. 


2 .  The  Error  Functional  and  Blending  System  of  Equations 

The  basis  of  the  FIB  analysis  is  the  minimization  of  an  error  functional, 
E  ,  over  the  two-dimensional  array  of  grid  points.  This  error  functional  for 
FIB-CLRX  is  given  as 
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In  the  error  functional,  the  <P  represent  the  blended  resultant  corrections 

u.n 

to  be  added  to  the  measured  normalized  radiance  anomalies,  T  .to 

v ,  n 

obtain  the  clear-column  radiance  components ,  3  .  The  entire  array  of 

"ft  ^ 

O y  ^  is  what  we  wish  to  arrive  at  by  blending  all  available  sources  of 
information. 


The  sources  of  information  are  combined  or  assembled  to  produce 

estimates  of  the  individual  corrections  <p  ,  u  and  L  N  .  Hence, 

U,n  v,r\  n  V,r\ 

estimates  of  the  individual  corrections  <p  are  formed  from  the  available 

U  ,n 

information.  The  more  valid  the  estimates  are  deemed  to  be,  the  higher  the 

square  of  the  disparity  between  the  assembled  values  and  the  blended 

resultant  values,  <p  ,  are  weighted  in  the  error  functional;  i.e. ,  the 

k  /  n 

higher  the  associated  weights  A  .  Similarly,  estimates  u  of  the 

v,n  '  v,n 

difference  between  the  corrections  for  two  adjacent  spots,  O  ,  -  <p 

v,nU  P,n 

are  formed.  Again,  the  square  of  the  disparity  between  the  assembled  values 
and  the  blended  resultant  values ,  <p*  ,  -  <p*  ,  are  weighted  as  to 

reliability  by  the  weight  ^  .  Note  that  since  we  cannot  spread  information 


across  the  lateral  boundaries  (i.e.,  spots  n=l  and  n  =  N),  the  weights  B 


v.O 


and  ^  are  always  set  equal  to  zero  where  they  occur  in  the  error  functional, 
and  in  the  blending  equations  to  be  derived  from  the  error  functional. 


Independent  estimates  of  the  entire  vector  of  corrections,  <P  ,  for 

~  n 

a  spot  may  be  made.  This  vector  may  be  written  in  the  normalized,  sign- 
adjusted  form  given  by 


N 

—  n 


IfSn1  'fin1 


(A2)* 


Note  that 

v 


N 


2 

v  ,n 


-  1  . 
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The  square  of  the  disparity  between  the  blended  resultant  value,  <P*  n  ,  and 

the  individual  component,  L  N  ,  of  the  total  adjustment  vector  are 

n  v ,  n 

weighted  as  to  reliability  by  the  weight  C  \V 

n  v ,  n 

It  is  important  to  note  the  kind  of  information  that  is  assimilated 

into  each  of  the  terms  in  the  error  functional.  Information  pertaining  only 

to  a  single  correction  goes  into  the  assembled  estimate  of  0 ^  ^  . 

Information  on  the  horizontal  gradients  or  spot-to-spot  differences  in  <P  for 

a  given  channel  goes  into  the  assembled  estimates  of  n  .  Information 

pertaining  to  the  whole  vector  of  corrections,  <P  n  ,  rather  than  just  the 

individual  corrections,  0  ,  goes  into  the  assembled  estimates  of  the 

V  ,r\ 

vector,  N  ,  and  of  the  magnitude,  L  . 

~  n  n 

The  blending  equations  are  obtained  by  setting 


9  E 


d<p* 

v,n 


0 


(A3) 


for  all  (f,n)  .  This  procedure  produces  the  general  equation 


+ 


AV,n  (^y.n  °t/,n)  Bp,n(<0i/,n+l  ^y.n  %,n) 

Bp,n-l(<Py,n  ~  0P,n-l  ~  %,n-l) 


+ 


C  W 
n 


(  O * 

V  ,n\  P,n 


L  N 
n 


v  ,n) 


0  . 


(A4) 


We  also  wish  to  make  that  adjustment  to  the  Ln  parameters  which 
will  minimize  the  error  functional.  By  setting 


3E 

BL 

n 


0 


(A5) 
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for  all  n  we  obtain 


Z 


N  W 
U,n 


( <p * 

v,n  \  v,  n 


L  N, 
n  v 


.») 


0 


(A  6) 


or 


L 

n 


H 


Z 

v 


N  W  (P 
v  ,n  v  ,n  v,r\ 


where 


H 

n 


W 


v,n 


(A7) 


(A8) 


Therefore,  we  may  write  the  blending  system  of  equations  (A4)  in  the 
equivalent  form 


v ,  n  ( 

<p* 
v  ,n 

<P  ) 

v  ,n/ 

v,n  \ 

* 

<p*  + 

y,ni-l 

„-i) 

v,n-l 

\  V,n 

-  <D* 

v,n-l 

-  U 

v, n-1 ) 

:  w 
n  v, 

( <p* 

,n V  v,n 

-  N  H 

u,n  n 

z 

to 

N  W  <P*  ) 

uj,n  to,n  u),n/ 

=  0 

• 

(A9)* 

★ 

The  same  form  for  the  blending  system  of  equations  is  obtained  by  substituting 
the  form  for  Ln  given  by  Eq.  (A7)  into  the  error  functional  Eq.  (Al)  and 
applying  Eq.  (A3). 
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This  system  may  be  expressed  in  the  stencil  format  given  in  Fig.  A1 .  The 
stencil  on  the  left-hand  side  of  Fig.  A1  shows  the  coefficients  that  multiply 
the  O  elements.  The  corresponding  forcing  function  elements  are  shown 
on  the  right-hand  side  of  Fig.  A1 . 


Blending  equations  in  stencil  format.  The  left -hand-side  shows  the 
coefficients  that  multiply  the  <p*  elements.  The  right-hand-side 
shows  the  corresponding  forcing  function. 


Since  the  blending  system  of  equations  is  a  set  of  simultaneous 
equations  which  are  linear  in  the  O*  elements,  they  may  be  expressed  in 
matrix  form  as 


It  will  be  propitious  to  use  the  following  ordering  for  the  unknown  (P ' 
elements: 


C*  i 

—  n-1 


(A  10) 


* 

O  = 


where 


(All) 


(A1 2) 


All  13  y  q  -  13^  ^  =  0,  as  noted  earlier.  The  coefficient  matrix  M  may 
then  be  written  in  partitioned  form  as 

M  = 


The  qeneral  submatrix  A  is  given  by 

;»n 


S>.n 

R, 

1 ,  n 

Rl 

1  » n 

R, 

1 ,  n 

Rl’.n 

Vn 

‘4, 

S0 

2  ,  n 

R2,„ 

R2 

2 ,  n 

R2,„ 

R2 ,  n 

R2 ,  n 

R3,n 

Rln 

S3,n 

r! 

3 ,  n 

R3,n 

3  ,n 

r73. 

A 

fa  n 

rJ 

4,n 

Rln 

R3 

4  ,n 

S4.n 

ir! 

4  ,n 

R? 

4 ,  n 

R,7 

4 ,  n 

R5,n 

R2 

5  ,n 

5 ,  n 

rJ 

5 ,  n 

S5,n 

Rr 

5 ,  n 

R2 

5  ,n 

R6 ,  n 

R2 

6,n 

4 

6,n 

Rf 

6 ,  n 

R6,n 

S6,n 

R7 

6 ,  n 

R7 

7,n 

R7 

7 ,  n 

R7 

7 ,  n 

R7 

7  ,n 

R7 

7,  n 

R7 

7 ,  n 

S7,n 

where 

A  + 

v,n 

B 

V  ,n 

+  B 

V ,  n 

-1  +  C 

W„ 

n  v ,  n 

fl  -H 
\  n 

N2 

i/,n 

and 

uu 

R 

v  ,n 

-  C 

n 

W„ 
v  ,n 

H  N 
n  V , 

N 

n  w,n 

W 

UJ,  n 

(A16) 


(A17) 


(A18) 


l)j  V 

Note  that  R.  =  R  so  that  A  is  symmetric. 
t/,n  yj,n  «n 

B  „  =  B  „T  =  0  . 
u,0  V,N 


Also,  here  again  all 
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The  general  .submatrix  B 

«n 


is  given  by 


B  = 
ss  n 


-B 


1  ,  n 
0 

0 

0 

0 


-B, 


2  ,  n 
0 

0 

0 


-B 


0 

3  ,n 
0 

0 


-B 


0 

0 

4,n 

0 

0 


-B 


0 

5 ,  n 

0  -B 

0 


0 


0 


6,n 


0 

0 


0  -B. 


7,n 


(A19) 


The  matrix  M  of  (A10)  may  therefore  be  described  as  a  symmetric  band 
% 

matrix . 
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3.  Methods  of  Solution  of  the  Blending  System 
3 . 1  Exact  Solution 

The  solution  to  the  blending  system  of  equations,  (A  1 0 ) ,  may  be 
written  as 


<P*  =  t  (A20) 

where  M  1  is  the  inverse  of  the  coefficient  matrix,  M  ,  and  is  a  full 
matrix,  in  general.  The  elements  of  M  will  be  given  by  M  ^  =  fm.  .  1  . 

»  RS  i ,  J 

A  linear  subscript,  p  ,  may  be  used  to  indicate  the  ordering  of  the  <p* 
elements  specified  by  (All)  and  (A12);  i.e.. 


<P 


* 

P 


(A21) 


where 


p  =  (n-l)Pm  +  v 
p  =  1 , 2 , . . . ,  Pm  , 

V  =  1,2, ,  (A22) 

and 


n  =  1 ,2  ,  . .  .  ,N 


From  (A20)  the  solution  for  the  arbitrary  element,  <p*  , 

P 


may  be  v/ritten  as 


Z 

i 


m  F 
PJ  J 
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(A23) 


where  the  subscript  j  ranges  over  all  values  of  p  from  1  to  p  ;  i.e.,  all 

n 

elements  of  the  forcing  function  contribute  to  all  elements  of  <p*  .  Trom 
(A14)  we  may  rewrite  (A23)  as 


(A24) 


The  object  of  the  blending  procedure  is  to  determine  not  only  the 

resultant  elements  of  cp  but  also  the  associated  resultant  weights,  A 

~  P,n 

A  basic  tenet  of  the  implicit  method  of  structure  blending  is  that  the 
resultant  blended  information  at  a  given  location  is  the  sum  of  the 
independent  information  assembled  at  that  location  plus  the  independent 
ambient  information  or  information  propagated  into  the  location  from 
surrounding  points  by  the  higher-order  parameters;  i.e. , 


A*  <P*  =  A  (p  +  A3  pa 

p  p  P  P  P  P 


(A25a) 


B* 

P 


B  M 
P  P 


Ba  M3 
P  P 


(A25b) 


A 

P 


(A2  6a) 


B 

P 


+ 


(A2  6b) 


The  first  term  on  the  right-hand  side  of  (A25)  and  (A2  6)  represents  assembled 

information  and  the  second  term  represents  ambient  information.  Eq.  (A24) 

may  be  rewritten  in  the  form 


( P 


* 

P 


m  A  <p  + 
P»P  P  P 
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(A27) 


where  X..*  refers  to  the  summation  in  (A24)  less  the  term  m  AO. 

p  p ,  p  p  p 

Since  X..*  contains  no  explicit  reference  to  A  O  ,  then  by  comparison 
v  p  P  P 

with  (A25a)  it  refers  to  the  independent  ambient  information  at  point  p  while 
the  first  term  refers  to  the  assembled  information.  Consequently,  we  may 
divide  both  sides  of  (A25a)  by  A^  and  equate  the  first  term  on  the  right-hand 
side  of  the  resulting  equation  to  that  of  (A27)  to  obtain 


A  0 

— ^ ^  =  m  AO 

A*  P,P  P  P 

P 


(A28) 


or 


(A29) 


•fa 

Hence,  the  blended  A  elements  are  given  by  the  reciprocal  of  the  diagonal 
elements  of  the  inverse  of  the  coefficient  matrix,  M  ,  of  (A10). 


The  blended  first  differences  may  be  written  as 


* 

P 


O 


p  +  V 


m 


0 


* 

P 


(A30) 


for  p  =  1 , 2 ,  . .  . ,  p  -  v  .  The  B  elements  may  be  obtained  by  employing 
m  m 

the  same  procedure  used  to  obtain  the  A  elements.  The  following  sequence 

of  equations  results: 
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M*  *  I 

P 


m  ,  ,,  |  A ,  O  .  -  B .  fX  .  +  B  fl  .  ) 

P  +  v  .  I  J  ]  i  }  j - p  j-v  / 

j  m,j  '  ■'mm' 


-  ?  mP.j(A)°>  -  Vi +  V*_  V-  ) 


m  m 


^  =  B  M  I  -  m  ,  +  m  ,  ,  +m  -m 

P  P  P  p  +  P  p  +  P  ,P  +  p  P,P  P 

1  m,p  mm 


(A31) 


,P  +  Pm 
m  J 


+  X  * 


(A3  2) 


K 


+  m  ,  -2m, 

p  p  +  p  ,p+P  p  + 

m  m 


■Vp] 


(A3  3) 


If  A  is  interpreted  as  the  reciprocal  of  the  unresolved  variance 
P 


in  (p  ,  then 
P 


<P 


(A3  4) 


Similarly  we  may  write 


2 

o  * 
p 


(A35) 


We  have  taken  advantage  of  the  fact  that  both  M  and  M  are  symmetric 

to  combine  the  m  .  and  m  .  terms. 

P,P+Pm  P-^  »P 
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Combining  (A29),  (A33),  (A34)  and  (A35)  we  arrive  at 


+  v 

m 


2 


m 

P 


i-U 

m 


P 


(A3  6) 


The  sum  of  the  variances  on  the  right-hand  side  of  (A3 6)  would  be  the 
variance  associated  with  the  first  difference  if  O  ^  and  <£>*+l,  were 


independent.  Hence,  the  subtraction  of  the  term  2  mp+li, 
the  interdependent  nature  of  the  lesultant  <P*  elements. 


m'P 


indicates 


The  solution  by  matrix  inversion  has  both  advantages  and  drawbacks. 
The  advantages  are  that  it  is  exact,  it  can  be  employed  for  any  form  of  the 
error  functional  which  produces  a  non-singular  coefficient  matrix  M  and  it 
yields  the  blended  weights  as  well  as  object  parameter  values.  The  main 
drawback  is  that  the  matrix  inversion  routines  may  be  prohibitively  time- 
consuming  compared  to  some  iterative  techniques. 


3 . 2  Point  Successive-Over-Rolaxation 

Point  successive-over-relaxation  (SOR)  could  be  employed  to 
iteratively  solve  (A10)  for  the  <P*  elements.  This  technique  is  well-known, 
and  well-documented.  Unlike  the  exact  matrix  inversion  method,  however, 
the  formulation  is  dependent  upon  the  form  of  the  error  functional  and 
resulting  stencil  for  the  blending  equations.  Solution  by  point  SOR  of  a 
blending  system  of  the  form  given  by  (A10)  and  a  perturbation  method  for 
approximating  the  A*  elements  are  discussed  in  f  1  ] . 


Holl,  Manfred  M.  and  Bruce  R.  Mendenhall,  1971;  "Fields  by  Information 
Blending,  Sea-Level  Pressure  Version",  Final  Report,  Project  M-167, 
Contract  No.  N66314-70-C-5226  (Fleet  Numerical  Weather  Central), 
Meteorology  International  Incorporated,  Monterey,  California,  66  pp. 
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3 . 3  Line  Successive-Over-Relaxation 


The  matrix  M  may  be  written  as 

ss 


M  =  A  +  B 

S3  S3  S3 


(A3  7) 


where 


A 

S3 


diag  A  ,  A 
«  1  «  z 


A.. 

S3  N 


(A38) 


and  B  is  given  by 


O 


b: 


ii 


o 

& 


rt 

I  2 
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O 

£3 


o 

BT 

Sn-i 


In-i 
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--  : - -  LTT" - 


The  blending  system  equations  may  then  be  rewritten  as 
(a  +  B^0*  =  F 


(A3  9) 


A  <P*  +  BT  .  (0*  _  f  B  (p* 
*»n  —  n  s»n-l  -n-1  «*n  — n+1 


n  =  1 ,2, . .. ,N 


where 


Equation  (A40)  may  be  solved  for  o  *  as 

—  n 


<P*  =  A  F  -  BT  .  O*  -  B  <p*  .1 

~n  «n  ~n  »  n-1  ~  n-1  «n  ~n+ll 


The  solution  by  line  relaxation  is  defined  by 


(e*„)(rtl>  *  (**„  (r’)+  [2  n  (r+1/z)  -  (2;)(r)' 

where  the  superscript  (r)  refers  to  pass  number  and  where 


(A40) 


(A41) 


(A42) 


*  \  (r+1/2) 


2*„) 


r  . 

v  <r)  , 

(r) 

(  An)~ 

1  I  F  -  BT  (0* 

\  (r+i)  _  /*  > 

1  (r+1)  j 

\«n/ 

1  — n  n-1  \  —  n-1 

/  »n\~-nHy 

'  J 

The  upper  superscript  (r)  applies  to  the  term  in  brackets  for  the  first  half 
of  each  pass  in  which  either  the  even  or  odd  numbered  n  are  solved  for. 
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*  -•  % 


The  lower  superscript  (r+1)  applies  to  the  term  in  brackets  for  the  second 
half  of  each  pass  in  which  the  remaining  (either  odd  or  even)  n  are  solved 
for.  The  method  converges  for  0  <  w  <  2  where  1  <  uu  <  2  implies 
overrelaxation . 

When  the  A  matrices  are  positive-definite  the  inversion  specified 
s»n 

by  Eq.  (A43)  may  be  obtained  explicitly  by  a  triangular  decomposition  of  the 

A  matrices  into  the  product  of  an  upper  and  a  lower  triangular  matrix, 

»  n 

followed  by  Gaussian  elimination.  If  a  large  number  of  passes  are  to  be 

used  and  the  weights  comprising  the  matrix  elements  do  not  change 

from  pass  to  pass,  then  a  significant  saving  in  computer  time  may  be 

[2  ] 

realized  by  this  decomposition. 

Note  that  solution  by  line  SOR  will  be  possible  as  long  as  the 
form  of  the  error  functional  yields  a  set  of  blending  equations  which  may 
be  put  into  the  form  given  by  (A10)  in  which  the  coefficient  matrix  M  may 
be  written  in  symmetric,  block  tri-diagonal  form. 


3.4  Method  Used  in  Development  Version  of  FIB-CLRX 

The  line  SOR  technique  has  been  employed  in  the  current  version  of 
FIB-CLRX.  After  the  explicit  solution  by  matrix  inversion  has  been  performed 
for  each  line  in  each  half  pass  of  the  iteration,  a  new  weight  Cn  is  computed 
using  the  relations  (A7)  and  (A50)  (to  be  discussed  in  Section  4).  The  new  C 

n 

is  used  in  the-  next  iteration  for  that  line.  Ten  passes  through  the  line  SOR 
technique  with  u)  =  1.2  give  adequate  convergence. 


r  2 1 

A  more  detailed  discussion  of  line  SOR  is  given  in  the  Appendix  to  "The 
Sea-Level  Pressure  Analysis  Capability,  FIB/SLP:  Improvements, 
Refinements,  and  Generalizations",  by  William  F.  Weigle  and  Bruce  R. 
Mendenhall,  Technical  Summary,  M-209,  Contract  No.  N66314-74-C-1528 
(Fleet  Numerical  V'eather  Central),  Meteorology  International  Incorporated, 
Monterey,  California,  November  1974,  37  pp.  plus  Appendix  A. 
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4. 


Assembly 


The  sole  source  of  information  for  this  developmental  version  of 

P 

FIB-CLRX  is  the  field  of  normalized  measured  radiance  anomalies,  H 

v,n 

from  a  scan  of  24  spots  of  7  channels  each.  From  these  anomalies, 

estimates  of  the  O  ,  M  ,  N  and  L  and  the  associated 
w,n  v,n  v,  n  n 

reliabilities  Ati  ,  B  ,  and  C  are  made.  Two  spots  have  been 

v ,  n  v ,  n  v ,  n  n 

omitted  on  each  end  of  the  scan  to  simplify  computational  loops.  Hence, 
for  blending  the  n  index  runs  from  3  to  22  and  the  v  index  from  1  to  7, 
producing  a  blending  system  of  140  simultaneous,  linear  equations.  The 
assembly  procedure  will  now  be  described. 

The  following  first-difference  vectors  are  formed: 


Ip 

~  n 

= 

~  n+1 

-  7  ‘ 

—  n 

n  = 

1,2, 

. . . ,  N-1 

(A44a) 

*2„ 

=  3P 
—  n+1 

-  SP 
—  n-1 

n  = 

2,3, 

.  . .  ,N-1 

(A44b) 

-  T p 

~n+2 

-  T P 
~  n-1 

n  = 

2,3, 

...,N-2 

(A44c) 

I 

From  each  of  these  first-difference  vectors  the  normalized  vectors 


The  N  vectors  surrounding  each  point  are  then  used  to  provide 

an  estimate  of  the  adjustment  vector  N  at  the  point  as  follows: 

—  n 


y  w 1  n  1 


N 


1/2 


(A46) 


i  ii^ii2  iAi  Ai 

where  W  =  u>  |  |  and  the  uj  ,  6  and  N  are  specified  in  Table  A1 

n  n  —  n 

below.  An  estimate  of  the  variance  associated  with  the  component  N  is 

u,n 

given  by 


2  w1  (nJ.  -  N  ) 

~  \  v  ,r\  v  ,n  / 


V,n 


(A47) 


W 


The  reliability  associated  with  the  component  N  for  blending  is  given 

”»n 

by 


W. 


V.n 


1 


a 


2 

v,n 


+ 


e 

v 


where  e  is  a  minimum  variance  associated  with  the  N 

v  v,n 


(A48) 
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Table  A1 


First  Differences  and  Weights  used 

to  Estimate  Adjustment  Vector  N 

~  n 


.  ^  i  ^i  Spots  surrounding  nth  spot  used 

1  ~n  ~~  n  to  estimate  difference  at  n _ 

n-2  n-1  n  n+1  n  (-2 


The  have  been  estimated  as  follows: 

€l  =  0.001 

e2  =  0.001 

c  3  =  0.001 

=  0.005 

4 

e  =  0.015 

J 

=  0.025 

b 

=  0.030  (A49) 

Note  that  for  a  given  N  ,  the  larger  L  is  the  larger  the  magnitude  of  the 

~n  n 

vector  difference  to*  -  L  N  is  and  the  less  the  associated  weight,  C  , 

~  n  n  ~  n  n 

should  be.  Hence,  the  weight  is  given  by 


C 


n 


(A50) 


where  the  minimum  variance  is  given  by  q  =  1 .  Hence,  the  combined  weight, 

C  W  ,  will  be  largest  where  there  is  little  disparity  between  the  first- 
n  v ,  n 

difference  estimates  used  to  produce  the  N  and  the  associated  magnitude, 

v,n 

L  ,  appears  to  be  small.  Conversely,  the  combined  weight  will  be  smallest 
n 

where  there  are  large  differences  between  the  first-difference  estimates  and 
the  associated  magnitude  appears  to  be  large. 

The  assembled  value  of  (0.  for  blending  is  given  by 

v,n 

PtJ  „  =  0  (A51) 

v,n 


for  all  v  and  n.  Consequently,  the  weight  A..  ,  should  reflect  how  close 

v,n 

to  the  clear  column  value  the  measured  radiance  anomaly  is.  A  method  of 


weighting  has  been  used  which  appears  to  give  satisfactory  results  for  all 
cases  except  those  where  a  uniform,  cloud  layer  exists  over  many  adjacent 
spots.  The  first  step  is  to  compute  at  each  point  a  measure  of  the  disparity 

P 

in  IT  between  adjacent  spots  according  to  the  formula 


°'2}2  *■  1  -  °'2}2 

°-4}2  *  |hi,n  -  I  -  0.2}2  (A52) 


for  V  -1 ....  ,7  and  n  -  3 ,  .  . . ,  22  .  The  square  of  negative  quantities  in  (A52) 
is  set  equal  to  zero.  A  tentative  weight  field  is  then  given  by 


(0.25  F  +  0.2)-1 
'  v,n  > 

(0.2  F  +  0.05  F  ,  +  0 . 2) _1 

'  v ,n  v-I , n  / 


v  =  1 

V  =  2,...,7  (A53) 


and  the  weight  field  for  blending  by 


'k  if  &  5  1 

v,n  v  ,n 


v.n 


0  if  A  <  1 
V,n 


(A54) 


The  first  difference,  u  ,  may  be  written  as 

—  n 


n  =  3, ... ,21 


n  =  1,2,22,23,24 


(A55) 
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*  p  * 

Since  S  =  S  +  P  is  what  we're  trying  to  determine  by  the  blending 
operation,  the  second  term  on  the  righi-hand-side  of  (A55)  is  initially  set 
equal  to  zero.  The  associated  weights  are  given  by 

n  =  3  ,  .  .  . , 21 

n  =  1 ,2,22,23,24  .  (A56) 


★ 

The  first-guess  field  of  p  lor  the  blending  operation  is  set  equal 

to  zero  for  all  (v,n)  where  A„  3  1.0.  Where  A  =0.0,  the  first¬ 
ly  n  v  ,n 

guess  value  is  the  difference  between  the  linearly  interpolated  value  of 

p 

S’  ,  using  the  nearest  adjacent  bracketing  spots  with  A  £  1.0,  and 

v.n  p  v,n 

the  reported  value  of  S 

v ,  n 


5. 


Recycling 


After  assembly  and  blending  the  <P*  and  A*  elements  reflect  the 
assimilation  of  all  available  information  according  to  th'  formulation  of 
the  error  functional  and  blending  system  of  equations.  Hence,  cp  *  and  A* 
may  be  used  to  compare  each  individual  piece  of  information  for  consistency 
with  all  other  information.  These  comparisons  will  then  be  the  basis  of  a 
reevaluation  of  the  weights  of  the  individual  information  elements  used 
during  assembly. 

Recycling,  then,  is  the  iterative  process  involving  the  above- 
mentioned  reevaluation,  assembly  and  blending  stages  of  the  FIB  methodology. 
It  begins  after  the  initial  assembly  and  blending  operations.  Its  purpose  is 
to  assess  more  accurately  the  independent  worth  of  the  individual  information 
elements  that  were  used  during  assembly. 

Recycling  could  continue  ad  infinitum;  however,  there  is  a  trade-off 
between  increasing  computer  time  and  decreasing  return.  In  previous  FIB 
applications,  three  cycles  (one  initial  assembly  and  blending  and  two 
reevaluation  cycles)  have  been  adequate.  Consequently,  we  have  employed 
three  cycles  in  FIB-CLRX  with  good  results  to  date.  The  details  will  now  be 
given. 

The  following  changes  are  made  to  the  weights  and  values  for  the 
second  and  third  cycle  assembly  and  blending: 


S® 


ten 


(A57) 
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i 


1.0 


A  (j)  A  (j) 

i  A  for  AVJ' 
v,n  V,  n 


0  for  A(j)  <  1.0 

V.n 


(A58) 


where  (  )  1  refers  to  cycle  number  2  or  3.  Note  that  this 

formulation  permits  information  which  was  rejected  during  a  previous 
cycle  to  be  reinstated. 


0  as  in  cycle  1 , 


(AS  9) 


25  . 


uyn  +  fK  for n  =  3,..., 21 


for  n  =  1,2,22,23,24 


(A60) 


where  f  =  0.8  , 


y  x*(i-d 
v  >n 


n  =  13 


12 

y  j*(j- 

4%  v,n 
n  -  3 


(A61) 


3P  + 
v,n  v,n 


The  term  fK^  in  (A60)  is  an  estimate  of  the  term  in  (A55)  that  was 
set  equal  to  zero  during  the  first  cycle.  It  is  inclusion  of  this 
term  that  permits  increasing  the  associated  weights  B 

v  ,n 
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The  calculation  of  the  N  by  (A46)  and  W  by  (A47)  involve 

.  —  n  v ,  n  A  . 

the  weight  factors  W1  and  the  individual  estimates  N1  .  The 

~  n 

reevaluation  consists  of  measuring  the  consistency  between  the 
a  i 

estimates  and  the  blended  adjustment  vector 


I  •  O 


N 

—  n 


-n 


< P 


* 

<P 

~n 


V 


-n 


and  using  this  measure  to  reevaluate  the  coefficient  |  |  ^ 

~n 

in  the  definition  of  the  weight  factors  W1.  Equations  (A46)  and 
(A47)  are  then  applied  with  the  new  weight  factors.  The  following 
formulas  are  used  to  reevaluate  the  coefficients: 


where  the  symbol  •  refers  to  the  dot  product.  When  the  initial 
estimates  and  the  blended  adjustment  vectors  are  parallel,  the 
dot  products  equal  one.  Where  they  are  not  parallel,  the  values 
range  between  one  and  zero.  Negative  values  are  set  equal  to 
zero  since  estimates  which  are  more  than  90°  out  of  phase  with 
the  blended  adjustment  vector  are  deemed  worthless. 
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The  net  result,  then,  is  that  in  both  (A46)  and  (A47) 

those  initial  estimates  of  N  which  are  more  consistent  with  the 

~  n 

blended  analysis  are  reinforced  during  reevaluation  while  those 
which  are  less  consistent  have  their  effect  reduced. 
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